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The  objective  of  this  experiment  was  to  show  the  relationship 
between  output  frequency  shifts  and  the  absorption  cell  magnetic  field  during 
microwave  power  fluctuations  on  a  commercial  (EG&G  Frequency  Products; 

Model  RFS-10)  rubidium  frequency  clock.  This  experiment  was  conducted  at  the 
Aerospace  Corp.  frequency  laboratory  in  El  Segundo,  Ca. 

Changes  in  the  clock  output  frequency  relative  to  the  10  MHz 
reference  were  monitored  as  changes  in  the  absorption  cell  magnetic  field  (0.01 
to  1.0  Gauss)  and  changes  in  the  cell's  interrogation  microwave  power  were 
varied  from  the  nominal  factory  specification  (P0)  at  four  discrete  values  (+1.3,  - 

1.1,  -2.2,  -4.0  dB).  The  data  presented  here  represents  the  difference  between 
two  1000-second  frequency  averages,  calculated  as  the  difference  in  output 
frequency  at  the  average  high  power  (fH)  and  that  at  the  average  low  power 
(fL)  and  normalized  to  the  nominal  output  (i.c.  (f^  -  fj_ )/l 0  MHz).  The  maximum 
normalized  output  frequency  shift  of  1.4X10"1®  occurred  at  -4dB.  The  maximum 
normalized  frequency  change  for  the  +1.3  and  -1.1  dB  data  is  about  3  x  10"11  and 

for  the  -2.2  dB  data  is  about  6  x  10'11.  As  a  function  of  the  change  in  microwave 

power,  the  maximum  frequency  change  is  about  2.6x10  1  1  per  dB  of  power 
change.  There  also  existed  several  magnetic  field  values  which  resulted  >n 

minimum  output  frequency  change.  This  occurred  at  absorption  cell  magnetic 

field  values  of  9  and  27  milli-Gauss,  where  the  nominal  setting  is  137  milli-Gauss. 

A  minor  design  change  in  the  electronics  package  is  recommended  to 
control  input  microwave  power  to  within  0.01  dB.  Use  of  low  cel!  magnetic  field 
values  to  control  frequency  shifts  cannot  be  made  at  this  tine  until  a  larger 

sample  of  clocks  are  tested. 


I.  / ntroduction 


Atomic  frequency  standards  are  very  stable  and  accurate  frequency 
sources.  Space  applications  rely  heavily  on  accurate  clocks  for  navigation  and 
communication  applications,  which  in  turn  are  based  on  frequency  standards. 
(20:537-539) 

This  section  provides  an  introduction  to  the  components  of  atomic 
standards,  to  the  measurement  approaches  used,  and  a  brief  description  of  what 
makes  atomic  standards  operate.  Also,  a  description  of  the  operating 

characteristics  of  the  specific  rubidium  frequency  standard  used  in  this 
experiment  will  be  provided.  At  the  conclusion  of  this  section  there  will  be  a 
discussion  of  why  the  test  results  are  of  importance  and  what  impact  they  could 
have  on  future  rubidium  dock  design  and  operation. 

Frequency  is  a  quantity  closely  related  to  time,  sometimes  referred  to 
as  its  inverse.  In  „ctualiiy,  frequency  is  derived  from  the  observations  of 
periodic  events,  which  occur  at  regular  intervals  (of  time).  If  constant,  this 
interval  is  called  the  period  (T)  and  the  frequency  of  the  oscillation  (D ) , 

measured  in  Hertz,  is  the  inverse  of  the  period: 


\>  =1/T 
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One  hertz  is  one  period  (oi  cycle)  per  second.  An  outstanding  feature  of 
frequency  and  time  is  the  high  accuracy  of  this  basic  definition  and  precision 
of  the  measurements,  which  during  the  last  twenty-five  years  have  progressed 
to  extremely  precise  levels.  (18:2) 
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The  basic  unit  of  time  is  the  second  (s).  Since  1967,  one  second  ts 
defined  as  "the  duration  of  9,192,631,770  periods  of  the  radiation  corresponding 
to  the  transition  between  the  two  hyperfine  levels  of  the  ground  state  of  the 
Cesium- 133  atom"  (14:4),  where  the  accuracy  of  this  period  is  good  through  the 
tenth  significant  digit.  Therefore,  the  frequency  of  the  Cesium  oscillations  are 
9.192631770  GHz.  Because  of  this  type  of  accuracy,  cesium  clocks  have  been 
widely  used  as  standards  for  precise  time,  as  well  as  used  for  test  equipment,  to 
evaluate  new  types  of  frequency  standards.  If  cesium  is  so  accurate,  then  why 
arc  other  types  of  standards  being  investigated?  Compared  to  rubidium  clocks, 
cesium  is  on  the  average  15  times  larger,  10  times  heavier  and  requires  twice 
the  power.  However,  cesium  still  out  performs  rubidium  in  terms  of  long-term 
accuracy  by  more  than  an  order  of  magnitude.  Because  of  these  reasons,  and 
the  need  to  conserve  space  and  power  for  satellite  operations,  increased 
research  into  improved  rubidium  clock  accuracy  has  been  accomplished  over 
the  past  several  years. 

Accuracy,  Reproducibility ,  Stability 


The  performance  of  frequency  standards  and  clocks  is  described  in 
terms  of  accuracy,  reproducibility  and  stability.  Accuracy  is  the  degree  of 

conformity  of  a  measured  and/or  calculated  value  to  a  specified  value. 
Reproducibility  is  the  ability  to  reproduce,  independently,  a  previous  frequency 
value.  Stability,  which  is  the  area  of  most  interest  in  frequency  standard 
development,  is  the  time  domain  behavior  (i.c.,  variations  in  frequency  over 
time)  of  a  process.  Measured  frequency  stability  depends  also  on  the  exact 
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measurement  procedure  used  to  conduct  an  experiment  or  evaluation. 

Frequency  stability  depends  on  the  tested  frequency  standard  and  the 
measurement  system  and  procedure  used.  The  frequency  stability,  or 
instability,  of  the  reference  measurement  can  only  be  ignored  if  it  is 
substantially  better  than  that  of  tb^  tested  device.  (14:5)  For  this  reason,  an  HP 
5061  Cesium  standard  was  used  as  the  reference  frequency  standard,  which 
exhibits  very  good  stability  characteristics  (10'14  Hz). 

In  the  time  domain,  the  acceptable  measure  of  stability  is  the  Allan,  or  two 
sample,  variance.  Frequency  stability  is  expressed  as  the  two-sample  (Allan) 
variance,  Oy2(t): 


Gy2(t)=  l/2[yk+l  -  yk]2  1‘: 

where  the  difference  in  the  brackets  denotes  an  average  over  a  number  of  data 
points  and  yk  is  the  average  fractional  frequency  over  time  of  the  data  point  k. 
(see  Appendix  A).  The  Allan  variance  is  used  since  frequency  measurements  in 
these  clock  standards  do  not  exhibit  a  normal  distribution.  The  quantities 
associated  with  frequency  measurements  are  the  number  of  measurements  (N), 
the  duration  of  each  measurement  (T ),  and  the  regular  measurement  interval 
(T).  The  variance  Oy2(x)  of  the  process  is  then  defined  in  terms  of  the  sample 
variance  Oy2(N,T,X)  as 


cry2(X)  =  lim  Cy2(N,T,X) 
N->infinity 
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when  the  limit  exists.  Even  after  elimination  of  any  initial  offset  and  linear 
drift,  it  has  been  found  that  the  sample  variance  Oy^fN.T.t)  depends  on  all  three 
variables  N,  T  and  t.  Therefore,  the  sample  variance  defined  in  Equation  1-2  has 
become  known  among  specialists  in  the  time  and  frequency  field  as  the  "Allan- 
variance".  (18:18) 

Resonators  (14:6-9) 


The  performance  of  a  frequency  standard  is  given  by  its  frequency 
determining  element.  It  determines  the  frequency  by  its  resonance  behavior. 
Some  examples  of  these  include:  a  rod  clamped  at  one  end  which  can  vibrate,  a 
block  of  solid  material  which  can  contract  and  expand  (thus  vibrate),  or  an 
antenna  (a  dipole)  where  the  distribution  of  electric  charges  can  oscillate  back 
and  forth.  These  devices  have  in  common  that  they  can  vibrate  or  oscillate  if 
they  are  excited,  by  mechanical  or  electrical  stimulation.  These  excited 
oscillations  will  cease  gradually  with  a  decay  time  which  is  determined  by  the 
losses  in  the  individual  resonators.  Determination  of  the  decay  time  (td)  of  the 
oscillations  will  describe  the  quality  of  the  resonator.  There  is  also  a  resonance 
frequency  at  which  all  resonators  oscillate  at  the  greatest  intensity.  On  both 
sides  of  the  resonance  frequency,  the  response  of  the  resonator  will  lessen  until 
it  ceases  to  respond.  There  exists  a  frequency  interval  around  this  resonance  in 
which  the  resonator  response  is  relatively  strong  and  is  called  the  resonance 
linewidth  (W).  The  relationship  between  linewidth  and  decay  lime  is  given  by: 

W=l/td  1-4 
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The  fractional  linewidth  is  determined  by  dividing  the  linewidth  W  by 
its  resonance  frequency  (t)0).  The  quality  factor,  Q,  is  more  commonly  used  as 

part  of  the  quality  specifications  for  frequency  standards  and  is  the  reciprocal 
of  the  fractional  linewidth: 


Q=Uo/W 
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This  is  graphically  represented  in  Figure  1.1.  Large  values  of  Q  are 
required  in  frequency  standards  to  achieve  accuracies  in  the  range  of  10'^  to 
10' I4,  which  implies  that  the  center  frequency  of  the  resonance  can  be 
determined  within  a  very  small  fraction  of  its  width. 


Figure  1.1  Decay  Time,  Linewidth.  and  Q-Value  of  Resonator  (14:7) 
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Quartz  Crystal  Resonator  The  quartz  crystal  is  a  mechanical  resonator 
which  needs  an  external  source  to  excite  its  oscillations.  This  occurs  due  to  the 
piezoelectric  properties  of  a  certain  class  of  crystals.  Compression  or  expansion 
of  the  crystal  generates  a  voltage  across  the  crystal  and  the  application  of  an 
external  voltage  across  the  crystal  causes  the  crystal  to  expand  or  contract 
depending  on  the  polarity  of  the  voltage.  (14:9) 

An  atomic  frequency  standard  is  developed  from  the  use  of  a  crystal 
resonator  as  the  external  voltage  source,  and  by  adding  an  electronic  amplifier 
and  feedback  system  back  into  the  atomic  resonator.  Its  output  frequency  is 
determined  by  the  quartz  crystal  resonator  whose  frequency  is  determined  by 
the  physical  dimensions  of  the  crystal  combined  with  the  crystalline  quartz 
properties.  The  nominal  output  frequencies  from  the  crystal  oscillator  are  5 
and/or  10  Mhz.  However,  electronic  multipliers  can  be  used  to  provide  any 
desired  output  frequency.  Therefore,  a  crystal  oscillator  with  a  required  output 
of  10  MHz  has  a  crystal  thickness  of  just  a  few  tenths  of  a  millimeter.  The  output 
of  the  crystal  oscillator  is  what  is  used  to  drive  other  systems  or  used  to  measure 
the  accuracy  of  the  standard  as  a  whole.  (2:3) 

Crystal  resonators  have  Q-values  which  are  typically  in  the  range  of 
Q=105  to  107.  However,  due  to  external  effects  such  as  temperature  and 
acceleration,  the  crystal  resonance  accuracy  deteriorates  more  rapidly  than  is 
required  for  space  applications.  If  they  are  calibrated  against  high  accuracy 
frequency  standards,  such  as  atomic  resonators,  they  maintain  this  calibration 
according  to  their  long-term  stability  performance.  (14:16-18) 

Atomic  Resonators  The  basis  of  commercially  available  atomic 
frequency  standards  are  resonance  in  atoms  at  microwave  frequencies,  usually 
in  the  range  of  1  to  100  GHz.  An  explanation  of  atomic  resonance  is  compared  to 
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a  magnetic  dipole  antenna.  The  dipole  antennae  can  be  separated  into  two  types, 
(1)  receiving  antennae,  which  absorb  energy  from  a  field  at  their  resonance 
frequency,  and  (2)  transmitting  antennae  which  radiate  energy  at  their 
resonance  frequency.  The  atoms  which  act  as  a  receiving  antenna  are 
identified  as  being  in  the  "lower  state"  and  those  acting  as  a  transmitting 
antenna  are  in  the  "upper  state".  This  is  represented  in  Figure  1.2.  (14:18) 


Receiver  Transmitter 

-  upper  state  - 


photon 


lower  state 


photon 


Absorption 


Emission 


Figure  1.2  Atomic  Resonance  as  Compared  to  Magnetic  Dipole  Antenna 


The  nominal  state  of  most  gases  have  the  total  number  of  upper  state 
atoms  closely  equal  to  the  number  of  lower  state  atoms.  In  actuality  though,  the 
probability  P(£)  of  finding  an  entity  in  one  of  its  variety  of  energy  states  £, 

for  a  system  containing  entities  at  the  same  thermal  equilibrium,  is: 

P(£)  =  A  e  'e/e0  1-6 

where  A  is  the  state  degeneracy  constant  and  £0  is  defined  as  the  quantity  kgT. 
The  value  kg  is  known  as  the  Boltzman  constant  equal  to  1.38  x  10' 16  erg/deg  and 
T  is  the  thermal  equilibrium  temperature  of  the  physical  entities  of  the  system. 
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system.  P(e)  is  referred  to  as  the  Boltzman  probability  distribution.  At  room 
temperature  (-300° K)  and  a  Ae  of  about  4.5x10' 17  ergs  between  the  Rubidium 
energy  states,  there  is  less  than  a  1%  probability  difference  between  finding 
atoms  in  the  upper  or  lower  energy  states.  (13:58-61) 

If  the  giS  is  placed  in  an  external  magnetic  field,  and  stimulated  by 
oscillations  at  the  atomic  resonance  frequency,  the  atoms  will  resonate.  Nearly 
half  of  the  atoms  absorb  energy  from  the  field  and  the  other  half  emit  energy  of 
an  equivalent  amount  to  the  field.  The  gas  as  a  whole  acts  as  if  it  has  no 
resonance,  although  the  individual  atoms  are  resonating.  In  order  to  observe 
the  atomic  resonance,  the  relative  amounts  of  each  of  the  states  of  the  atoms 
must  be  changed.  The  upper  or  lower  state  has  to  be  in  the  majority, 
corresponding  to  a  net  absorption  or  emission  of  energy.  This  process 
determines  the  design  of  an  atomic  resonance  device.  The  frequency  of  the 
output  is  associated  with  the  transitions  between  these  two  energy 
levels.  (14:19) 

State  Selection  There  are  two  methods  to  accomplish  the 
change  in  the  relative  number  of  the  two  kinds  of  atoms:  1)  spatial  state 
selection  and  2)  optical  state  selection.  The  discussion  below  is  depicted  in 
Figures  1.3  and  1.4,  respectively. 

Spatial  state  selection  relies  on  the  actual  sorting  of  the  two  atomic 
states  into  different  directions  in  space.  One  of  these  states  is  then  used  as  the 
resonance  standard  and  the  other  is  discarded.  Usually  in  this  case  an  oven  is 
used  to  produce  an  atomic  gas  by  heating  it  and  allowing  the  atoms  to  exit  the 
oven  through  a  small  hole  which  forms  an  atomic  beam  in  an  adjacent  vacuum 
chamber.  This  atomic  beam  then  passes  through  a  magnet  which  causes  the 
separation  of  the  beam  into  two  separate  beams,  each  of  them  containing  only 
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one  of  the  atomic  states.  This  occurs  because  the  magnet  will  exert  a  force  or. 
these  atoms  which  will  act  in  opposite  directions  for  the  two  atomic  states.  This 
method  is  usually  used  in  the  design  of  Cesium  frequency  standards.  (14:19-20) 
The  optical  state  selection  technique  takes  advantage  of  the  fact  that 
the  atoms  have  more  than  one  resonance.  These  resonances  typically 
correspond  to  infrared  or  visible  light  frequencies.  One  of  these  resonances 
can  be  excited  by  shining  an  intense  light  of  corresponding  frequency  on  the 
atoms.  These  visible  light  wavelengths  for  the  ^^Rb  and  ^Rb  atoms  will  be 
discussed  in  Chapter  2  and  how  these  atoms  partially  overlapping  wavelengths 
are  used  to  take  advantage  of  a  phenomenon  referred  to  as  optical  pumping.  The 

optical  state  selection  method  is  usually  used  in  the  design  of  Rubidium 
frequency  standards  (14:20). 

Atomic  Interrogation  (14:20-21)  A  microwave  cavity  is  used 
to  observe  the  atomic  resonance  in  both  of  the  state  selection  techniques.  These 

cavities  have  microwave  resonances  which  are  determined  by  their  physical 
size.  The  electrical  losses  of  these  cavities  are  determined  by  the  electrical 
conductivity  of  the  cavity  material.  These  cavities  are  normally  shaped  as  a 
cylinder  or  rectangle,  depending  on  the  desired  effect.  In  order  to  observe  the 

atomic  resonance,  the  state  selected  atoms  are  placed  inside  the  cavity  and 
subjected  to  a  microwave  signal  corresponding  to  their  resonance  frequency. 
The  microwave  signal  will  change  the  the  relative  number  of  atoms  in  the  two 
states.  If  all  the  atoms  were  initially  in  the  upper  state,  then  after  the 
microwave  signal  has  acted  upon  them  there  will  be  lower  state  atoms  produced. 
If  the  frequency  of  the  external  microwave  signal  is  exactly  at  the  atomic 
resonance  frequency,  this  transfer  of  atoms  from  one  state  to  the  other  reaches 
a  maximum.  The  center  of  the  atomic  resonance  is  then  found  by  monitoring 
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the  number  of  atoms  in  one  of  the  two  states  and  by  varying  the  microwave 
frequency  until  the  number  of  atoms  in  the  lower  state  reaches  a  maximum,  or 
the  number  of  atoms  in  the  upper  state  a  minimum.  A  proportional  electric 
signal  is  then  derived  which  is  fed  back  to  the  oscillator  generating  the 
microwave  signal.  This  automatic  servo  mechanism  keeps  a  crystal  oscillator 
locked  to  the  atomic  resonance. 

Signal  Detection  (14:22)  The  detection  of  the  effect  of  the 
microwave  signal  on  the  atoms  is  accomplished  in  three  different  manners;  1) 
atom,  2)  optical  and  3)  microwave  detection.  The  optical  detection  method  is  the 
most  commonly  used  in  present  rubidium  standards  and  will  be  the  only  method 
discussed  here.  The  other  methods  of  signal  detection  are  discussed  in  detail  in 
Reference  14. 

If  the  atoms  are  optically  state  selected,  a  photo-detector  is  placed  in  a 
position  such  that  the  light  which  has  interacted  with  the  atoms  is  detected. 
Since  the  light,  which  enters  the  cavity  at  the  same  wavelength,  removes  atoms 
from  one  of  these  states,  the  light  intensity  will  change  when  it  exits  the  cavity. 
This  occurs  since  the  atoms  in  the  cavity  will  absorb  the  energy  of  the  light 
photons  and  decrease  its  intensity  as  the  light  strikes  the  photodctector.  The 
light  at  the  photo-detector  is  then  a  measure  of  the  number  of  atoms  in  one  of 
the  states,  due  to  the  light  intensity  change.  Figure  1.2  graphically  depicts  how 
the  energy  of  the  light  photon  can  move  atoms  from  one  state  to  another.  As 
was  seen  in  Figure  1.3,  the  optical  wavelengths  of  the  two  popular  types  of 
rubidium  have  a  very  close  visible  spectrum  and  this  natural  phenomenon  is 
taken  advantage  of  in  the  design  of  rubidium  clocks  and  will  be  explained  in 
Chapter  II  in  the  optical  pumping  section. 
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Fractional  Frequency  (23:4-5)  In  providing  information  regarding 

frequency  accuracy  or  stability  for  a  given  device,  the  most  common  method 

has  been  termed  fractional  frequency.  Fractional  frequency  is  a  dimensionless 
quantity  which  determines  an  error  from  some  known  or  standard  frequency 
value.  In  most  applications  of  frequency  standards,  the  output  frequency  used 
to  drive  other  devices  is  not  the  atomic  clock  resonance  frequency.  Electronic 
frequency  multiplication  or  division  is  applied  to  provide  a  more  standardized 

value,  such  as  5  or  10  Mhz.  Because  of  these  possible  applications,  absolute 
frequency  values  are  rarely  used.  The  term  fractional  frequency  (y)  has  been 

adopted  to  provide  a  common  basis  of  world-wide  discussion,  and  is  given  by  the 
relationship 


y  =  (Um-'V/'>n  1-7 

where  \)m  is  the  measured  frequency  value  and  \)n  is  the  nominal  value 
expected.  The  term  "fractional",  though  in  common  use  throughout  the  English 
speaking  world,  is  being  discouraged  by  the  Committe  Consultatif  Internationale 
de  Radio  (CC1R),  because  of  the  difficulties  in  translating  that  term  into  certain 
other  languages.  The  CCIR  prefers  the  use  of  "relative”  or  "normalized".  Results 
of  the  experiment  conducted  here  will  be  provided  in  terms  of  fractional 
frequency,  but  the  above  comment  should  be  noted. 

Rubidium  Clock  Principle  of  Operation 

The  rubidium  system  is  composed  essentially  of  an  optical  package 
acting  as  a  frequency  reference,  and  of  a  frequency-locked  loop  for  locking  the 
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frequency  of  a  crystal  oscillator  to  the  optical  package  reference  as  depicted  in 
Figure  1.5.  The  optical  package  consists  of  a  ^Rb  iamp(  a  8 5 hyperfine  filter, 
a  ^Rb  absorption  cell,  a  microwave  cavity  and  a  photodetector.  A  solenoid  is 
used  to  create  a  weak  magnetic  field  inside  the  absorption  cell,  which  is  usually 
encased  in  magnetic  shielding  to  attenuate  the  effect  of  external  magnetic  field 
fluctuations.  (2:1259) 

Light  from  the  ^Rb  iamp,  filtered  by  the  ^Rb  hyperfine  filter,  is 
focused  on  the  absorption  cell  containing  87Rb  and  the  buffer  gas.  The  lamp 
consists  of  a  small  glass  envelope  approximately  1  cm  in  diameter,  containing 
about  0.5  mg  of  ^Rb  and  argon  at  a  pressure  of  approximately  2  Torr.  It  is 
excited  by  an  oscillator  delivering  a  power  of  1  Watt  at  a  frequency  of  80  MHz 
and  heated  to  a  temperature  of  125°  C.  The  spectrum  of  light  emitted  consists 
mostly  of  the  D1  and  D2  lines  of  ^Rb  as  seen  jn  Figure  1.5.  As  a  fact  of  nature, 
one  of  the  ^Rb  wavclength  lines  is  very  closely  coincident  with  one  of  the 
87Rb  wavelength  lines.  (2:1260) 

LAMP  HYPERFINE  MAGNETIC  MICROWAVE 


Figure  1.5  Rubidium  Frequency  Standard  Block  Diagram  (2:1260) 


1.13 


The  emitted  light  is  allowed  to  pass  through  the  hyperfine  filter 
containing  the  8^Rb  p]us  the  buffer  gas.  The  closely  spaced  optical  spectrum 
lines  of  the  two  gases  is  why  rubidium  is  used  in  atomic  frequency  standards, 
and  is  exploited  for  filtering  the  87Rb  spectrum.  The  filter  cel!  contains  an 
argon  buffer  gas  which  slightly  shifts  the  8^Rb  spectrum  towards  the  red 
wavelength  and  broadens  the  lines  which  make  up  the  overlap  to  upper 
wavelength.  This  increases  the  efficiency  of  the  filter.  When  the  resulting 
light,  after  going  through  the  filter,  enters  the  absorption  cell  containing 
additional  87 Rb,  only  atoms  in  the  ground  state,  which  are  resonant  with  the 
state  of  interest,  are  then  affected.  These  atoms  are  excited  to  a  higher  energy 
state.  (2:1261) 

Through  spontaneous  emission  or  relaxation  by  collisions  with  the 

buffer  gas  in  the  absorption  cell,  the  atoms  fall  back  to  the  ground  state.  Since 
the  process  of  excitation  is  continuously  done  from  the  lower  ground  level  tFl), 

the  effect  is  to  increase  the  population  of  the  level  F2  at  the  expense  of  the  FI 

level,  which  in  turn  provides  an  atomic  population  imbalance  between  these 

two  ..tates.  This  effect  is  known  as  optical  pumping.  (9:599-600) 

Referring  to  Figure  1.6,  the  net  result  of  the  optical  pumping  from 
the  lamp  is  to  overpopulate  F2  with  respect  to  FI.  Then  a  microwave  radiation  at 

the  frequency  corresponding  to  the  F2  to  FI  (mf=0  to  mf=0  )  transition  is  applied 

(theoretically  6. <*346826128(5)  GHz),  moving  the  atoms  back,  toward  equilibrium. 
When  the  microwave  frequency  is  exactly  on  the  resonance  value  of  a  specific 

absorption  cell,  then  the  total  light  transmitted  through  the  absorption  cell  to 

the  photodctcclor  is  at  a  minimum.  (2:1263) 

A  constant  magnetic  field  is  applied  to  the  absorption  cell  to  keep 
these  ground-state  energy  levels  apart.  It  was  determined  by  Robert  H.  Dicke  of 
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Princeton  University,  that  if  the  container  held  only  rubidium  atoms,  the  Rb 
irregular  thermal  motion  would  cause  a  condition  known  as  a  Doppler  shJt. 

This  effect  shifts  the  resonant  frequency  by  a  different  amount  for  each  atom, 
thus  smearing  out  the  resonance  line.  (5:78) 

The  accuracy  of  the  standard  depends  on  the  relationship  between  the 
magnetic  field  and  the  position  of  the  levels.  Frequency  stability  of  the  standard 

depends  on  the  size  of  the  magnetic  field  perturbations  and  the  ability  to  control 
the  field.  In  order  to  excite  the  transitions,  a  constant  magnetic  field  must  exist 
at  the  site  of  the  atom  and  must  be  parallel  to  the  RF  field  which  excites  the 
transitions.  (2:1260) 

EG&G  Model  RFS-/0  Description  (19:1-1  to  4-8) 

The  RFS-10  rubidium  clock  manufactured  by  EG&G  Frequency 
Products  in  Salem,  MA.  is  basically  a  crystal  oscillator  locked  to  an  atomic 
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Figure  1.6  Stale  Transitions 


1.15 


reference.  A  cross  sectional  view  is  shown  in  Figure  1.7.  The  rubidium  physics 
package  serves  as  a  passive  discriminator,  producing  an  error  signal  that  varies 
in  magnitude  as  a  function  of  the  difference  in  frequency  between  the  applied 
rf  excitation  and  the  atomic  resonance.  The  rubidium  frequency  is  about  6835 
MHz,  with  a  Q  of  10^.  The  error  signal  is  a  result  of  an  FM  audio  frequency  (-500 
Hz)  applied  to  the  microwave  excitation,  which  causes  variations  in  the 
transmission  of  light  through  the  resonance  cell.  A  photodetector  senses  this 
response,  which  has  a  typical  signal-to-noise  ratio  of  70  dB,  in  a  1  Hz  bandwidth. 

The  physics  package  operation  is  supported  by  an  exciter  for  the 
rubidium  lamp  and  temperature  controllers  for  the  lamp  and  cell  ovens.  The 

error  signal  is  processed  by  a  servo-amplifier,  which  generates  a  voltage  that 
controls  the  frequency  of  a  crystal  oscillator.  This  oscillator  produces  the 
output  and  also,  via  a  synthesizer  and  multiplier  chain,  drives  the  physics 
package.  The  overall  scheme  is  that  of  a  frequency  locked  loop. 

The  RFS-10  resonance  cell  consists  of  two  devices,  a  filter  cell  and  an 
absorption  cell.  This  arrangement  generally  gives  better  performance  than  if 
the  two  functions  were  performed  by  a  single  integrated  cell. 

The  RFS-10  rubidium  physics  package  is  contained  inside  a  1.3  inch 
diameter  cylindrical  magnetic  shield  about  2.6  inches  long. 

Dimensions  (19.1-4) 

The  RFS-10  Rubidium  Frequency  Standard  is  a  miniaturized  high- 
performance  unit  measuring  2-3/4  inches  by  2-3/4  inches  by  4  inches.  It 
features  low  operating  power  (12  watts  at  25  degrees  C),  fast  warm-up  (5  minutes 
at  25  degrees  C),  and  ruggedized  design  for  military  applications. 


ouitu  cover  l  motherboard 


i 


T 


Figure  1.7  RFS-10  Cross  Sectional  View 
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The  RFS-10  also  features  double-shielded  construction  for  low  magnetic 
sensitivity.  The  lamp  and  cell  ovens  have  low  thermal  mass  to  aid  in  fast  warm¬ 
up. 

The  RFS-10  has  two  ovens  that  control  the  temperature  of  elements  inside 
the  rubidium  physics  package.  One  oven  heats  the  rubidium  lamp  and  the  other 
heats  a  microwave  cavity  which  contains  the  Filter  and  absorption  resonance 
cells  and  also  holds  a  crystal  oscillator. 


Physics  Package  Description  (19:4-3) 

The  physics  package  assembly  is  comprised  of  the  following  major 
elements:  1)  lamp  exciter,  2)  lamp  oven,  3)  microwave  cavity,  4) 
photodetector,  5)  crystal  oscillator  and  6)  preamplifier.  A  cross  sectional  view  of 
the  physics  package  is  depicted  in  Figure  1.8. 

The  lamp  oven  contains  the  rubidium  lamp.  The  microwave  cavity 
contains  the  filter  and  absorption  cells  and  a  microwave  multiplier.  The  heart 
of  the  RFS-10  unit  is  the  physics  package  described  in  the  following  pages. 

The  physics  package  includes  two  lenses  in  the  optical  path.  The  first, 
located  in  front  of  the  absorption  cell,  collimates  the  lamp  output  to  give  a 
uniform  light  intensity  distribution.  The  second  lens,  located  behind  the 
absorption  cell,  focuses  the  light  on  the  photodetector  to  give  good  collection 
efficiency.  The  photodetector  is  a  custom  EG&G  device  having  low  noise. 

The  rubidium  lamp  is  excited  with  an  rf  power  of  about  0.5  watt  at  105 
MHz.  This  energy  is  applied  to  the  lamp  by  operating  it  inside  the  field  of  an  rf 
coil  which  excites  optical  emissions  from  an  ionized  plasma.  The  lamp  provides 
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the  Rb  resonance  radiation  required  for  optical  pumping. 

The  step-recovery  diode  (SRD)  microwave  multiplier  is  incorporated  into 
the  temperature-controlled  cavity.  Its  design  is  both  simple  and  efficient  since 
it  works  as  a  suaight  x76  multiplier,  it  is  driven  by  a  +20  dBm  signal  at  about  90 
MHz  and  produces  about  100  micro- watts  at  the  rubidium  resonance  frequency 
of  approximately  6835  MHz. 

The  VCXO  board  is  located  inside  the  physics  package  assembly.  It  consists 
of  a  20  MHz  voltage-controlled  crystal  oscillator  (VCXO)  that  is  locked  by  the 
servo-amplifier  to  the  rubidium  resonance  frequency  and  produces  the  output 
signal  from  the  RFS-10.  The  VCXO  uses  a  20  MHz  fundamental  mode  SC-cut 
crystal.  This  crystal  type  is  chosen  for  its  small  size,  fast  warm-up,  and  low 
vibration  sensitivity.  It  is  mounted  on  the  1.25  inch  diameter  VCXO  board  that  is 
attached  to  the  photodetector  end  of  the  microwave  cavity  and  thereby  shares 
the  80  degrees  C  stabilized  thermal  environment  of  the  cavity  oven. 


Problem  Statement  and  Experiment  Approach 

There  is  a  close  bond  between  the  interaction  of  the  atomic  magnetic 
moment,  the  external  C-field  to  which  it  is  exposed,  and  the  power  associated 
with  the  electromagnetic  wave  which  causes  the  hyperfine  transitions. 

As  evident  in  experiments  conducted  by  Andrea  DiMarchi  (12:54-58),  the 
cesium  frequency  standard  can  be  improved  in  long-term  frequency  stability 
by  the  correlation  of  absorption  cell  magnetic  field  values  versus  the 
microwave  power  used  to  excite  the  cesium  hyperfine  transitions.  His 
experiment  and  analysis  concluded  that  there  were  certain  magnetic  field 
values  which  minimized  any  frequency  shifts  in  the  output  of  the  cesium  clock 
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due  to  any  changes  in  microwave  power.  Variations  in  microwave  power  could 
result  in  external  effects  on  the  clock  electronics  which  would  make  it  difficult 
to  control  when  used  in  space  applications.  His  recommendation  was  to  adjust, 
or  tune,  the  c^>um  cavitv  to  one  of  several  specific  magnetic  field  values, 
which  in  turn  improved  the  long  term  accuracy  of  the  clock  (11:517-521). 

Since  similar  experiments  have  not  been  performed  on  rubidium  clocks, 
this  paper  will  parallel  the  experiment  conducted  by  DiMarchi,  except  using  a 
commercially  available  rubidium  standard,  provided  by  EG&G  Frequency 
Products.  Results,  analysis  and  conclusions  will  be  provided  as  to  the  outcome  of 
this  experiment. 

The  objective  of  this  experiment  is  to  gather  and  analyze  applicable  data 
which  will  provide  the  necessary  information  to  validate  theoretical  models  and 
answer  the  primary  question:  "Are  there  magnetic  field  values  which  improve 
(a  specific)  rubidium  clock  performance  when  the  microwave  power  is 
varied?" 

The  following  questions  will  be  investigated: 

1)  Does  the  theoretical  model  for  rubidium  hyperfine  transitions  match 

the  experimental  results?  What  is  the  error  associated  with  experimental  data? 

2)  The  average  magnetic  field  value  can  be  calculated  using  existing 

models.  Can  this  model  be  used  to  accurately  determine  the  magnetic  field  inside 
the  absorption  cell? 

3)  Is  there  a  linear  or  non-linear  relationship  between  Zeeman 

frequencies  and  coil  magnetic  field  changes?  Explain  the  relationship. 

4)  Are  there  magnetic  field  settings  within  the  operating  parameters  of 

this  rubidium  clock  which  minimize  output  frequency  variations  to  microwave 
power  variations? 
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5)  Is  there  a  relationship  between  changes  in  microwave  power  and 
changes  in  the  values  of  the  Zeeman  frequencies? 

6)  Is  there  a  relationship  between  changes  in  microwave  power  and  the 
shape  of  the  resonance  spike? 
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II.  Theory 


A  basic  understanding  of  the  theory  of  atomic  structure  and  the  stability 
of  atoms  will  be  discussed  to  explain  the  operational  phenomena  within 
Rubidium  frequency  clocks.  The  chemical  properties  of  all  elements  in  the 
periodic  table  give  rise  to  their  suitability  for  use  in  accurate  long  term  time 
measurement  applications.  Hydrogen,  consisting  of  one  single  electron 
orbiting  its  nucleus,  has  been  studied  in  detail  and  gives  rise  to  the  most 
accurate  timing  standards  available.  Other  elements  exhibit  a  similar 
characteristic  to  that  of  hydrogen.  This  similarity  is  the  presence  of  an 
unpaired  electron  orbiting  the  nucleus  of  the  atom,  and  are  known  as  the  alkali 
elements  in  the  periodic  table.  These  elements  have  one  single  electron 
orbiting  around  a  closed  subshell  of  paired  electrons.  The  resonance  of  the 
single  orbiting  electron  is  one  of  the  reasons  that  rubidium  clocks  are  used 
today.  (2:7) 

There  exists  two  naturally  formed  rubidium  atoms;  one  has  an  atomic 
weight  of  85  and  the  other  87.  Both  of  these  atoms  are  used  together  in  the 
design  of  rubidium  clocks.  The  standard  nomenclature  to  differentiate  between 
these  two  atoms  are  the  symbols  ^^Rb  and  87Rb.  In  the  case  of  87Rb,  the  nucleus 
is  composed  of  37  protons  (known  as  the  atomic  number),  and  50  neutrons.  It 
consists  of  five  shells,  or  layers,  with  a  total  of  37  electrons.  Each  electron  shell 
is  filled  through  the  determination  of  the  individual  quantum  numbers  n,l,  and 
m.  The  principal  quantum  number,  n,  will  determine  the  allowable  values  for 
the  azimimuthal  and  magnetic  quantum  numbers  1  and  m.  The  quantum 
numbers  can  take  the  following  values: 
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n  =  1,  2,  3 .  infinity 

1  =  0,  1,  2 . (n-1) 

m  =  -I,  -(1-1) .  +(1-1),  +! 

The  particular  conditions,  regarding  the  minimum  and  maximum  values  that 
these  quantum  numbers  can  take,  are  set  by  the  eigenvalue  equations  resulting 
from  Equation  2-2  when  the  solution  for  Equation  2-4  is  replaced  in  it.  (8:  53) 
The  representation  for  the  electron  distribution  for  87Rb  is  given 

through  the  following  series  arrangement:  1S^  2S^  2P®  3S^  3P®  3D1®  4S^  4P6 
5S1.  The  numbers  in  front  of  the  letters  indicate  which  orbiting  shell  the 
electrons  are  found  in  (related  to  the  principle  quantum  number  n);  the  letters 
represent  the  electrons  angular  momentum  (related  to  the  value  of  1)  ,and  the 
superscript  identifies  the  number  of  electrons  in  that  subshell  composition. 

Each  of  these  states  correspond  to  a  different  energy  level. 

The  87 Rb  electron  distribution  above  indicates  that  the  last  shell  (5S1) 
contains  an  unpaired  electron  as  seen  in  Figure  2.1.  Movement  of  this  electron 
through  some  type  of  stimulation,  allows  it  to  move  to  another  energy  state.  As 
will  be  seen  in  this  chapter,  the  movement  between  these  energy  states 
corresponds  to  different  frequency  values.  Figure  2.1  also  represents  an 

exaggerated  view  of  this  electron  and  the  next  two  energy  states  (5P)  in  which 
this  electron  could  be  moved.  The  subscripts  in  this  figure  represents  the 
designation  for  the  total  angular  momentum  due  to  electron-nucleus 

interactions,  and  will  be  discussed  in  more  detail.  The  D1  and  D2  designation  in 
the  figure  represents  the  optical  wavelength  associated  with  the  movement  of 
the  electron  between  these  two  stales. 

The  interaction  between  the  electron  and  the  nucleus  plus  the  effects 
which  occur  due  to  outside  stimulation  will  be  briefly  discussed  in  this  chapter 
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to  provide  a  better  understanding  of  what  is  occurring  at  the  atomic  level  to 
provide  the  results  of  this  study. 


D2=780. 


Figure  2.1  87Rb  Unpaired  Electron  Configuration 


Atomic  Spectra 

The  field  of  atomic  frequency  standards  is  characterized  by  a  level  of 
accuracy  that  is  not  seen  in  other  fields  of  physics  and  engineering.  This 
accuracy  relies  partly  on  the  existence  of  calculations  of  various  perturbations 
which  influence  the  properties  of  the  atoms  of  interest.  As  errors  in  the 
neighborhood  of  10' 14  are  being  investigated,  precise  measurements, 
procedures  and  calculations  need  to  be  performed.  A  brief  discussion  on  atomic 
spectra,  as  well  as  some  details  on  rubidium,  will  follow  to  orient  the  reader  as  to 
why  some  phenomena  occur  when  atoms  are  subjected  to  certain  stimuli. 
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Principal  Quantum  Number  An  electron  in  an  atom  moves  about  the 
nucleus  under  the  influence  of  the  Coulomb  attraction  between  the  electron  and 
the  nucleus.  The  potential  energy  of  this  closed  system  is  given  by: 

V(r)  =  -e2/  4  7C  £0  r  2-1 

The  Schroedinger  time-independent  wave  equation  is  determined  from  the 
following  relationship:  (10:7-4) 

((  -  h2  /  2  [l  V2)  -  (e2  /  4  it  e0  r»  u(r,  0,  <J>)  =  E  u(  r,  0,  <J>)  2-2 

where  u  (r,  0,  <J>)  is  the  wavefunction  in  spherical  coordinates,  independent  of 
time,  H  is  the  reduced  mass  expressed  as: 

H  =  Mp  m  /  (Mp  +  m)  2-3 

and  Mp  is  the  mass  of  the  proton  and  m  is  the  mass  of  the  electron.  In  expression 

2-2,  E  is  the  total  energy  of  the  atom  assumed  at  rest,  r  is  the  electron-proton 
distance,  £0  is  the  permittivity  of  free  space,  e  is  the  electron  charge  and  V2  is 

the  Laplacian  operator. 

The  solution  of  this  equation  can  be  obtained  by  expressing  V2  in  terms  of 
spherical  coordinates.  The  variables  are  separated  and  the  wavefunction 
u(  r,  0,0)  can  be  written  in  terms  of  Laguerre  polynomials,  R,  and  spherical 
harmonics,  Y.  Several  functions,  u(  r,  0,0),  characterized  by  the  quantum 
numbers  n,  1.  and  m  are  obtained: 
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u„lm<  r.  0.  4»  =  R„,(r)  Ylm(0,  $) 
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The  normalized  spherical  harmonic  Yim(0,0)  can  be  written  in  terms  of 
the  normalized  associated  Legendre  polynomial  P\m  as 

Ylm(6,  <j>)  =  [l/(27t)  >'2]  Plm  e  *m0  2-5 

The  principal  quantum  number,  n,  is  generated  by  the  solution  of  the 
radial  part  of  the  equation.  It  gives  rise  to  the  energy  levels:  (10:7-4) 

En  =  -  (^  e4>/(32  7t2  E02  h1  n2)  2-6 

Since  all  items  are  constants  in  the  energy  equation  except  for  n,  this  can 
be  reduced  to: 


En  =  -  Ej/n^  2-7 

where  Ej  is  the  atom's  ionization  energy.  For  Rubidium,  Ej  =  4.177  eV  (electron- 
volts)  where  1  eV  =  1.5956  x  10"  ^  ergs.  The  lowest  energy  (the  most  negative)  is 
that  in  which  the  principle  quantum  number  n  equals  1,  which  is  called  the 
ground  state.  (25:924-925) 

Each  energy  level  corresponds  to  a  stationary  state  in  which  the  atom  can 
exist  without  radiating.  Stationary  states  above  the  ground  state  (n  =  2,  3,  4,  ..,) 
are  called  the  excited  states.  An  electron  in  one  of  these  states  tends  to  make  a 
transition  to  one  of  the  lower  stationary  states.  (25:925) 


2.5 


An  electron,  initially  moving  in  an  excited  state  of  total  energy  Ej 
changes  its  motion  such  that  it  moves  into  a  lower  state  of  total  energy  Ef.  In 
the  transition,  the  atom  losses  energy  (Ej  -  Ef),  where  a  single  photon  having 
energy  hv  is  created  and  emitted  by  the  atom.  The  frequency  of  the  emitted 
radiation  \)  is  equal  to  the  difference  between  these  energy  levels  divided  by 
Planck's  constant,  or: 


\)  =  (Ei  -  Ef)/h  2-8 

Therefore,  \)  would  be  considered  a  constant,  only  dependent  on  t'.it  two  levels 
chosen  in  the  energy  scheme  of  the  atom.  (15:773) 

Equation  2-8  assumes  that  the  atoms  are  at  rest.  However,  atoms  are 
continuously  in  motion  and  their  emission  and  absorption  frequencies  are 
functions  of  their  velocity  relative  to  the  observer.  Since  the  atomic  velocities 
arc  random,  with  a  Maxwellian  distribution  (2:451),  the  net  effect  is  that  each 
atom  emits  or  absorbs  energy  at  slightly  different  frequencies,  referred  to  as 
Doppler  broadening.  The  net  frequency  result,  U0,  called  the  resonance 

frequency,  is  spread  over  a  broad  range.  The  half  power  frequency  width  of 
this  doppler  effect  is  given  by: 

A\)D  =  2  \)0  [  (  2  kB  T/ M  c2  )  In  2]>/2  2-9 

where  M  is  the  rest  mass  of  the  atom,  and  the  other  variables  have  been 
previously  defined.  In  the  case  of  optical  spectra,  ih„  frequency  width  Adb  can 
be  very  large  due  to  the  fact  that  U0  is  high.  In  the  case  of  the  ground-state 
hyperfine  transition  in  Rubidium  at  room  temperature,  the  width  is 
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approximately  3  KHz.  The  principle  behind  an  atomic  frequency  clock  is  to 
reduce  this  broad  range  of  frequencies  to  a  very  stable  value  with  extremely 
small  perturbations,  in  the  short  or  long  term,  from  the  resonance  value. 
Accurate  controlling  of  this  frequency  provides  the  timing  accuracy  and 
stability  required  in  modem  time  devices.  (13:468-472) 

The  following  sections  are  only  approximate  descriptions  of  the  classical 
quantum  mechanics  associated  with  a  multiple  electron  atomic  configuration, 
and  are  provided  to  aid  the  reader  in  visualizing  the  relationships  of  angular 
momentum,  magnetic  moments,  and  internal  atomic  magnetic  field  interactions. 
The  unpaired  electron  in  rubidium  is  the  one  of  interest  as  compared  to  the 
single  hydrogen  atom,  but  there  other  relativistic  effects  from  the  inner  shell 
electrons  which  contribute  to  its  behavior. 

Orbital  Angular  Momentum  (25:930-931)  The  motion  of  the  electron 
around  the  nucleus  in  these  varictv  of  orbits  possesses  an  orbital  angular 
momentum.  In  the  ground  state,  where  the  electron  orbit  is  circular,  the 
orbital  angular  momentum  is  zero.  The  magnitude  of  the  electron's  angular 
momentum,  L,  is  represented  in  terms  of  the  orbital  angular-momentum 
quantum  number,  /,  and  is  given  by: 

L  =  [/(/+  1)  ]1/2  h  2-10 

where  h  is  Planck's  constant  (h)  divided  by  2k.  The  value  of  /  may 
take  on  integral  values  starting  with  zero  and  continuing  up  to  (n-1),  where  n  is 
the  principle  quantum  number.  Each  value  of  /  corresponds  to  a  s'^te.  The  case 
of  l  equal  to  zero  is  where  the  electron  has  no  orbital  angular  momentum 
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relative  to  the  nucleus  and  is  considered  to  be  in  the  S  state.  The  next  values  of  / 


(1,2, 3,4...)  correspond  to  the  states  P,D,F,G .  respectively. 

When  a  single  electron  undergoes  a  transition  from  one  allowed  state  to 
another,  and  a  photon  is  emitted,  that  state  can  only  change  such  that  the  value 
of  /  changes  by  the  integer  1,  or  A/  equals  +/-  1. 

Orbital  Magnetic  Moment  (13:327-329)  Associated  with  the  orbital  angular 
momentum  vector  of  the  last  section  is  the  orbital  magnetic  moment  p/ .  From 

electromagnetic  theory,  the  electron  charge  circulating  in  a  loop  will  generate 
a  current.  In  addition,  such  a  current  loop  produces  a  magnetic  field  which  is 
the  same  at  large  distances  as  that  of  a  magnetic  dipole  located  at  the  center  of 
the  loop  and  oriented  perpendicular  to  its  plane. 

For  a  current  i  in  a  loop  of  area  A,  the  magnitude  of  the  magnetic  moment 
of  the  equivalent  dipole  is: 


M-/  =  iA 
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Figure  2.2  is  a  representation  of  the  angular  magnetic  moment  vectors. 
The  magnetic  moment  is  a  measure  of  the  strength  and  orientation  of  the  dipole. 
The  abbreviations  np  and  sp  relate  to  the  strength  of  the  positive  pole  (np)  and 
the  negative  pole  (sp). 

For  a  current  loop  produced  by  a  negatively  charged  particle,  the 
direction  of  H/  is  opposite  the  direction  of  the  angular  momentum  vector  L.  Also, 

the  ratio  of  the  magnetic  moment  to  the  magnitude  of  the  angular  momentum  is 
a  constant  known  as: 


M-/  /  L  =  g i  [Lft/h 
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Figure  2.2  Vector  Representation  of  Orbital  Angular  Momentum 
Vector  L  and  Associated  Magnetic  Moment  Vector  [Ll  (13:328) 

where  pg,  called  the  Bohr  magneton,  forms  a  natural  unit  for  the  measurement 
of  atomic  magnetic  moments.  The  quantity  g/  is  called  the  orbital  g  factor,  with 
a  value  approximately  equal  to  one. 

The  magnetic  moment  can  be  written  in  terms  of  the  orbital  angular 
momentum: 


H/  =  *  (g /  Hb  L  2-13 

The  ratio  of  of  these  two  vectors  does  not  depend  on  the  size  or  shape 
(circular  versus  elliptical)  of  the  orbit  nor  on  the  orbital  frequency. 

Equation  2-10  provides  the  magnitude  of  the  angular  momentum  as  a 
scalar,  but  the  direction  was  not  identified.  The  orbital  magnetic  quantum 
number  m/  yields  the  component  of  the  electron's  orbital  angular  momentum. 
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relative  to  some  arbitrary  direction  in  free  space,  which  will  be  called  z.  The 
values  of  m  /  range  from  +  /  to  -  /.  Figure  2.3  represents  the  magnitude, 

direction,  and  z-component  of  the  angular  momentum  vector  associated  with  an 
electron  in  the  D  state  (/  =  2). 

Since  the  value  of  the  z-component  of  the  angular  momentum  vector  can 
be  represented  in  terms  of  the  quantum  number  m /,  then  by  Equation  2-13  the 

magnetic  moment  vector  in  the  z-direction  can  also  be  represented  in  terms  of 
the  same  quantum  number: 

4/z  =  '  Si  M-B  ml  2-14 

As  seen  in  Figure  2.3,  Lz  =  nt[h  with  m  [  an  integer,  and  since  the 
direction  of  the  z  axis  was  chosen  arbitrarily,  then  the  quantum  mechanics 

ml  k 

2  }  2  h 

Z 

1  lh 

0  Oh 

-1  -lh 

-2  -2  h 

Figure  2.3  Representation  of  Orbital  Angular  Momentum  Vector  L 
and  Orbital  Magnetic  Quantum  Number  m/ for  D  State  (25:932) 
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requires  the  component  in  any  direction  of  the  orbital  angular  momentum  to  be 
an  integral  multiple  of  h.  This  is  ^nown  as  spatial  quantization. 

Atomic  Magnetic  Field  Effects.  (13:330-331)  If  the  dipole  represented  in 
Figure  2.2  is  subjected  to  an  external  magnetic  field  H  in  the  z  direction  ,  then 
there  will  be  a  force  acting  on  the  positive  and  negative  poles  (np  and  sp).  In 
this  external  magnetic  field  a  force  of  magnitude  spH  will  act  on  the  positive 
pole  (np).  with  an  equivalent,  but  oppositely  directed,  force  acting  on  the 
negative  pole  (sp).  There  is  no  net  translational  force,  but  there  will  be  a  torque 
applied  to  the  center  of  the  structure.  The  magnitude  of  the  torque  will  also 
depend  upon  the  angle  between  the  magnetic  moment  |X/and  the  vector  Hz. 

The  magnitude  and  direction  of  this  torque  can  be  expressed  through  the  cross 
product  of  the  two  vectors: 


T  =  \it  x  H 
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Therefore,  with  this  orientation  torque,  the  atomic  system  will  have  an 
orientation  potential  energy,  AE,  dependent  only  on  the  angle  between  the  two 
vectors.  This  AE  can  be  written  in  terms  of  the  dot  product: 


AE  =  -  •  H 
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This  orientational  potential  energy  is  minimized  when  the  two  vectors  are 
aligned.  Since  there  is  no  internal  mechanism  to  allow  for  a  dissipation  of 
energy,  then  the  magnetic  moment  will  prccess  around  the  external  magnetic 
field  vector  in  such  a  way  that  the  ar  1c  between  the  vectors  and  AE  will  remain 

constant. 
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This  precessional  motion  is  a  consequence  of  the  torque  always  being 
perpendicular  to  the  angular  momentum  vector  L.  Equation  2-16  can  be  written 
in  terms  of  L,  from  what  was  derived  in  Equation  2-13.  Therefore  L  and  | i./  will 
precess  about  H  with  an  angular  frequency  of  magnitude  COl  and  can  be 
represented  as: 


“L  =  M-B/  h)  H 
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This  is  known  as  the  Larmor  frequency  and  is  graphically  represented  in 
Figure  2.4. 

Electron  Spin  (13:334-338)  In  addition  to  the  angular  momentum  due  to 
an  electron  orbiting  the  nucleus,  it  also  has  an  intrinsic  angular  momentum  due 


Figure  2.4  Magnetic  Moment  Vector  \il  Precessing 
About  External  Magnetic  Field  H  (13:332) 
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to  the  existence  of  the  electron  spinning.  The  magnitude  of  this  spin  angular 
momentum,  S,  is  related  to  the  spin  quantum  numbers  s.  The  equation  for  the 
angular  momentum  magnitude  is  the  same  as  that  for  the  orbital  angular 
momentum  vector  (reference  Equation  2-10)  and  given  by: 

S  =  [  s  (s+1)  ]1/2  h  2-18 

The  relationship  between  the  spin  angular  momentum  and  the  spin 
magnetic  moment  is  the  same  as  before  and  can  be  written  as: 

4s  =  '  4B  /  h)  s  2-19 

where  gs  in  this  case  is  called  the  spin  g  factor,  which  is  approximately  equal  to 
two.  The  actual  value  of  gs  for  Rubidium  is  2.002331.  (2:24) 

Since  the  value  of  the  z-component  of  the  angular  momentum  vector  can 
be  represented  in  terms  of  the  quantum  number  ms ,  as  in  Equation  2-14,  then 

the  magnetic  moment  vector  in  the  z  direction  can  also  be  represented  in  terms 
of  this  spin  quantum  number: 


4sz  -  *  Ss  4b  ms  2-20 

Due  to  the  fact  that  the  axis  of  rotation  does  not  change,  it  can  be  seen 
that  the  z-component  of  the  magnetic  moment  for  a  single  electron  atom  will 
only  assume  two  values,  which  are  equal  in  magnitude  but  opposite  in  sign.  The 
spin  quantum  number  ms  can  only  have  the  values  +/-  1/2  and  s  can  only  have 

a  single  value  equal  to  1/2. 
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The  spin  angular  momentum  and  the  magnetic  moment  are  opposite  in 
direction,  as  were  the  orbital  angular  momentum  and  the  orbital  magnetic 
moment  (seen  in  Figure  2.4). 

The  atomic  system  will  also  have  an  orientational  potential  energy 
contribution  due  to  the  spin  magnetic  moment.  This  AE  can  be  written  in  terms 
of  the  dot  product: 


AE  =  -  •  H  2-21 

Total  Electronic  Angular-Momentum  (13:346-353)  The  orbital  angular 
momentum  L  plus  the  spin  angular  momentum  S,  form  the  total  electronic 
angular  momentum  J: 


J  =  (L+S)  2-22 

For  simplification  then,  the  quantum  number  mj  was  formed  to  represent 
the  sum  of  the  two  other  quantum  numbers,  and  ms,  which  will  be  reduced 

when  substituting  in  the  values  for  the  orbital  and  spin  angular  momentum 
vectors  into  Equation  2-22.  Since  the  maximum  possible  value  of  m[  is  /,  and  the 
maximum  possible  value  of  ms  is  s=  1/2,  then  the  maximum  possible  value  of  mj  is 

(/  +  1/2).  In  the  ground  state,  where  L=0,  it  can  been  seen  that  J  =  S,  therefore, 
mj  equals  +/-  1/2,  and  for  j  it  is  equal  to  1/2. 

The  electron  spin-orbit  interaction  energy  contributions  were  directly 
related  back  to  the  quantum  numbers  n  and  j  to  produce  the  equation  which 
identifies  the  atomic  fine  structure  for  hydrogen  atoms  only: 

E  =  -  (n  Z2  e4  /  2n2  h 2)  [1  +  (a2  Z2  /  n)(  1/0+1/2)  -  3/(4n))]  2-23 
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where  a  is  called  the  fine  structure  constant,  equivalent  to  approximately  1/137 
or  more  exactly  (7.29720  x  10-3)  (10:7-27). 

This  equation  determines  the  energy  levels  in  terms  of  the  electron  orbit 
and  received  consistent  results  with  Bohr's  findings,  when  the  electron’s  orbits 
were  assumed  to  be  circular  in  the  calculation  for  the  state  energy  values.  As 
can  be  seen  in  Figure  2.5,  the  number  of  energy  splits  at  each  level  is  related  to 
the  principal  quantum  number  n.  This  also  identifies  the  only  permitted 
energy  level  transitions  between  these  energy  splits,  which  is  consistent  with 
the  explanation  of  the  quantum  number  l  discussed  earlier.  These  fine  lines  are 
greatly  exaggerated  for  representation  purposes.  The  difference  between  these 
lines  at  each  energy  state  level  is  on  the  order  of  a 2.  The  arrows  indicate 

transitions  between  the  various  energy  states  allowed  which  produce  the  lines 
of  the  atomic  spectrum. 
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Figure  2.5  Splitting  of  the  n=2  Levels  by  the  Spin-Orbit  Coupling 

Plus  Relativistic  Effect 
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Zeeman  Effect(  13:443-448) 


It  has  also  been  observed  that  placing  an  atom,  such  as  the  alkali  atoms, 
in  an  external  magnetic  field  contributes  to  what  is  now  called  the  Zeeman 
effect.  This  effect  is  an  observation  that  the  spectral  lines  produced  due  to  the 
fine  structure  are  also  separated  into  several  components.  For  external  fields 
up  to  a  few  thousand  gauss,  this  splitting  is  proportional  to  the  strength  of  the 
magnetic  field,  but  is  small  in  comparison  to  the  fine  structure  splitting  of  the 
spectral  lines  of  the  atoms.  This  also  means  that  in  the  presence  of  a  magnetic 
field,  the  energy  levels  of  the  atom  which  contribute  to  the  spectral  lines  must 
then  also  be  split  into  several  components. 

Earlier  it  was  discussed  (page  2.12)  that  the  summing  of  the  electrons 
orbit  and  spin  angular  momentum  resulted  in  the  total  angular  momentum 
vector  J.  So  it  follows  that  in  the  presence  of  the  same  magnetic  field  H,  that  the 
energy  difference  can  also  be  represented  as: 


AE  =  -  IiH  H 
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where  [J.jj  is  the  component  of  the  total  magnetic  moment  along  the  direction  of 
H. 

The  Larmor  precession  that  existed  with  the  angular  momentum  and 
magnetic  moment  vectors  in  Figure  2.2  also  holds  true  for  the  spin  angular  and 
magnetic  moment  vectors.  In  this  case,  however,  L  and  S  precess  about  J 
because  of  the  precession  of  S  in  the  atomic  magnetic  field  associated  with  L, 
and  have  a  processional  frequency  w^.  This  precessional  frequency  is 

proportional  to  the  strength  of  the  atomic  magnetic  field.  The  total  magnetic 
moment  vector  H,  will  also  precess  about  J  with  the  same  precessional 
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frequency.  In  ihe  presence  of  an  external  magnetic  field  H,  there  will  also  be  a 
precession  of  p  about  the  direction  of  this  field  with  a  precessional  frequency 

proportional  to  the  strength  of  H.  If  the  external  field  is  weak  compared  to  the 
atomic  magnetic  field,  then  the  precession  of  ji.  about  H  will  be  slow  compared  to 

the  precession  about  J.  This  is  graphically  represented  in  Figure  2.6. 

All  cases  discussed  in  this  paper  assume  that  the  external  magnetic  field  is 
considerably  smaller  than  the  atomic  magnetic  field.  The  value  of  the  atomic 
magnetic  field  has  been  measured  to  be  on  the  order  of  10^  gauss. 

Since  ji  precesses  much  more  rapidly  about  J  than  about  H,  can  be 
evaluated  by  looking  at  |lj.  Therefore,  can  be  represented  in  terms  of  pB  and 
the  quantum  number  my.  The  resultant  average  energy  difference  can  then  be 
written  as: 


AE  =  pg  H  g  ntj 
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Figure  2.6  Electron  S  and  L  Angular  Momentum  Vectors 
Processing  about  J  which  Processes  about  Z 
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The  quantity  g  is  called  the  Lande  g  factor  and  is  composed  of  the  quantum 
numbers  s,  /,  and  j,  and  evaluated  through  the  expression: 


g  =1+1  j(j+l )  +  J(.r+1)  - 1(1+1)]  /  2 /  (j+ 1 )  2-26 

The  net  effect  of  this  is  graphically  represented  in  Figure  2.7.  The  vertical  lines 
indicate  the  only  allowed  state  transitions  between  energy  levels  and  follow  the 
earlier  defined  optical  selection  rule  (page  2.12)  of  Amy-0  and  +/-1.  The  S  state 

in  Figure  2.7  refers  to  the  ground  state  where  the  quantum  number  n=l, 
therefore  1=0,  and  the  subscript  1/2  refers  to  the  only  value  that  j  can  assume, 
which  is  1/2.  For  the  P  state,  n=l,  therefore  1  can  have  values  of  both  0  and  1. 
From  this,  the  subscript  of  1/2  and  3/2  relates  to  the  two  possible  values  of  j. 


No  magnetic  Magnetic  field 

field 


Figure  2.7  Zeeman  Splitting  of  Typical  S  and  P  Energy  Levels  (8:237) 
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Hyperfine  Structure  (13:449-452) 


In  addition  to  the  splitting  of  the  spectral  lines  associated  with  the  fine 
structure,  there  are  additional  separation  components  that  arise  due  to  the  effect 
of  the  nucleus  spinning.  These  lines  are  very  closely  spaced  and  are  on  the 
order  of  three  orders  of  magnitude  smaller  than  the  separation  characteristics 
of  the  fine  structure  of  spectral  lines.  These  small  components  are  referred  to 

as  the  hyperfine  splitting  of  the  energy  levels  of  the  atom. 

This  hyperfine  splitting  is  a  result  of  the  orientational  potential  energy 
of  the  magnetic  dipole  moment  due  to  the  atomic  nucleus  interacting  with  the 
magnetic  field  associated  with  the  motion  of  the  electrons.  The  nuclear 
magnetic  dipole  moment  is  due  t0  the  nuclear  sPin  an8ular  momentum,  I, 
which  is  derived  in  the  same  manner  as  the  electron  spin  angular  momentum 
vector,  S.  A  representation  of  both  the  electron  and  nucleus  angular 
momentum  vectors  is  depicted  in  Figure  2.8. 


Figure  2.8  Model  of  Electron  and  Nucleus 
Angular  Momentum  Vectors  (8:230) 
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The  magnitude  and  Z  component  of  I  are  represented  using  a  new 
quantum  number  associated  with  nucleus  spin  m,-,  and  is  derived  in  the  same 

manner  as  the  electron  spin  angular  momentum: 

l  =  [  <  0+1)  ]1/2  A  2-27 

therefore 

lz  =  m/  h  2-28 

The  equation  relating  the  nuclear  magnetic  dipole  moment  to  the  nuclear 
spin  is  given  by: 


Hi  =  (gf  h)  1  2-29 

such  that 

Hiz  =  giHNmi  2-30 

The  factor  g,-  is  also  on  the  order  of  1  (.995141  for  ^Rb)  (2:25),  but  is  positive  for 

some  nuclei  and  negative  for  others,  depending  on  the  spin  direction.  The 
quantity  |l(q,  which  is  called  the  nucleus  magneton,  is  a  factor  of  1836  smaller 
than  that  of  pg  (the  Bohr  magneton).  Since  the  nuclear  magneton  is  over  three 
orders  of  magnitude  smaller  than  the  Bohr  magneton,  it  follows  that  the  nuclear 
magnetic  moment  is  smaller  than  the  total  electron  magnetic  moment  by  the 
same  factor.  Consequently  the  orientational  potential  energy  of  the  nuclear 
magnetic  moment  (hyperfine  structure)  is  three  orders  of  magnitude  smaller 
than  the  orientational  potential  energy  of  the  atomic  magnetic  moment  (fine 
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structure),  since  both  magnetic  moments  experience  comparable  magnetic 
fields. 

The  interaction  of  the  nuclear  magnetic  dipole  moment  with  atomic 
magnetic  field  couples  the  nuclear  spin  angular  momentum  vector  I  to  the  total 
angular  momentum  vector  J  of  the  electrons.  Because  of  this  coupling,  a  new 
grand  total  angular  momentum  vector  is  generated,  F,  composed  of  the  sum  of  I 
and  J: 


F  =  I  +  J 
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Bo.ii  the  I  and  J  vectors  will  now  be  seen  preccssing  about  the  grand  total 
angular  momentum  vector.  Consequently,  when  the  hyperfine  structure 
interaction  is  evaluated,  m,  and  m j  are  not  used,  but  a  new  quantum  number,  ny, 

is  introduced.  The  relationship  between  the  magnitude  and  z  direction  is  the 
same  as  the  earlier  sums  seen  in  angular  momentum  and  grand  total  angular 
momentum.  Grand  total  angular  momentum  can  now  be  represented  as: 

F  =  [/(/+  1)  )1/2  h  2-32 


and  the  z  component  by: 


Fz  =  ny  /i  2-33 

Rubidium  Zeeman  Splitting 

The  effects  from  the  Zeeman  splitting  for  the  Rubidium  atom  discloses 
that,  in  the  visible  light  spectrum,  this  atom  produces  two  primary  wavelength 
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spectrums,  referred  to  as  the  D1  and  D2  lines.  They  are  represented  by  taking  a 
closer  look  at  the  frequencies  associated  with  the  transitions  between  the  Pj  /2 
and  P3/2  to  the  ground  state  Sj/2-  Figure  2.3  can  be  expanded  and  represented 
by  Figures  2.9  and  2.10.  These  figures  represent  the  state  transition  frequencies 
which  sum  together  to  give  the  two  optical  wavelengths  represented  by  the  D1 
and  D2  lines  (XI  =  794.76  nm  and  X2  =  780.02  nm).  (2:1260) 

The  same  evaluation  can  be  conducted  for  the  8^Rb  atom  to  determine  its 
optical  wavelengths.  As  can  be  seen  in  Figure  2.11,  one  of  the  wavelengths  for 
*^Rb  partially  overlaps  one  of  the  wavelengths  for  87Rb.  This  close  overlap  at 
the  higher  wavelength  provides  a  phenomenon  in  the  operation  of  atomic 
clocks  known  as  optical  pumping.  The  term  optical  pumping  refers  to  the 
redistribution  of  atoms  among  their  fine  or  hyperfine-structure  levels  by 
means  of  light.  As  part  of  the  design  of  rubidium  clocks,  a  Filter  cell  is  usually 
found  in  between  the  absorption  ceil  and  the  rubidium  lamp.  This  filter  cell  is 
filled  with  85Rb  atoms  and  an  argon  buffer  gas.  Because  the  D1  lines  do  not 
overlap  totally  for  85 Rb  and  87 Rb,  a  buffer  gas  is  used  in  the  filter  cell.  This 
buffer  gas  (usually  argon  or  nitrogen)  under  pressure  can  shift  or  broaden  the 
85Rb  lines  to  provide  for  a  good  overlap  of  the  87Rb  D1  line.  (5:73) 

Using  argon  as  the  buffer  gas,  the  following  approximate  values  for  the 
85 Rb  shift  and  broadening  parameters  have  been  reported:  A  shift  of 
-18  MHz/Torr  for  Dl,  -14  MHz/Torr  for  D2,  and  a  broadening  of  approximately 
10-20  MHz/Torr.  (2:1260) 

The  resultant  effect  of  overlapping  the  D2  lines  is  such  that  the  87Rb 
visible  light,  as  it  passes  through  the  filler  cell,  has  the  energy  at  the  Dl 
wavelength  absorbed  by  the  now  shifted  8^Rb  Dl  wavelength.  As  the  light  exits 
from  the  filter  cell,  it  is  composed  mostly  of  the  energy  associated  with  the  D2 
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Figure  2.9  Zeeman  Transitions  Producing  87Rb  D2  Line 


Figure  2.1 1  Rubidium  Optical  Emission 

line.  As  it  enters  the  absorption  cell,  which  contains  the  87Rb  atoms,  the  light  is 
absorbed  from  only  the  one  unfiltered  wavelength  line,  D2.  Thus,  an  imbalance 
in  the  fine  or  hyperfine-structure  is  achieved.  (2:1262) 

As  can  be  seen  in  Figure  2.11,  the  primary  resonance  frequency  is 
approximately  6.835  GHz  for  87Rb  and  approximately  3.036  GHz  for  8^Rb.  Both 
atomic  structures  can  be  used  in  the  operation  of  a  Rubidium  atomic  clock. 
However,  because  of  its  higher  frequency,  87Rb  provides  better  timing 
resolution. 

Rubidium  Hyperfine  Structure  (2:1265-1268) 

It  was  seen  that  in  the  presence  of  an  external  magnetic  field  the  energy 
levels  experience  a  Zeeman  splitting  effect.  In  the  ground  state  (^Sj^),  it  has 

been  seen  that  a  splitting  occurs  into  two  sublevel  ground  states  due  to  the 
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nuclear  spin  effect  within  the  atomic  magnetic  field  produced  by  the  electron. 
These  two  states,  named  FI  and  F2,  were  determined  through  Equations 
2-32  and  2-33.  The  Zeeman  splitting  of  each  of  these  levels  is  due  to  subjecting 
the  atomic  structure  to  an  external  magnetic  field.  As  the  external  magnetic 
field  (B)  increases,  the  distance  between  these  hyperfine  levels  increase.  An 
expanded  representation  of  these  hyperfine  levels  is  seen  in  Figure  2.12. 
Following  the  state  selection  rule  defined  before,  transitions  between  these 
hyperfine  states  can  only  occur  between  levels  which  differ  by  0  or  +/-  1. 


mf=2 


mf=l 


mf=0 


mf=-l 


mf=-2 


mf=-l 


mf=0 


mf=l 


Figure  2.12  Hyperfine  Splitting  of  87Rb  Ground  State 
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In  reality,  however,  the  Zeeman  splitting  does  not  occur  linearly.  As  the 
B  field  increases  to  large  values  (>100  gauss),  the  hyperfine  lines  become  non¬ 
linear.  The  regions  used  in  rubidium  clocks  are  linear,  since  the  B  fields  are 
controlled  under  1  Gauss.  The  bold  line  between  the  F2;  mf=0  and  FI;  mf=0  states 

occur  at  a  frequency  equivalent  to  6.835...  MHz.  These  transitions  are  the  ones 
used  in  the  rubidium  atomic  clocks.  A  more  accurate  representation  of  these 
hyperfine  energy  levels  is  presented  in  the  next  section. 

There  are  eight  different  possible  transitions  that  can  occur  between  the 
FI  and  F2  ground  states.  However,  there  are  two  pairs  of  transitions  that  are 
equivalent  in  frequency  on  both  sides  of  the  F2  (mf=0)  to  FI  (mf=0)  transition, 

which  sum  to  produce  a  slightly  stronger  signal.  So,  exciting  each  of  these 
transitions  at  the  appropriate  frequency  and  plotting  these  transition  lines, 
only  seven  frequencies  (also  called  the  Zeeman  frequencies)  will  be  seen.  These 
Zeeman  frequencies  are  referred  to  as  K  or  <J  transitions  based  on  the  value  of 

Amf.  The  K  transitions  refer  to  those  frequencies  associated  with  a  Amf=0,  and  O 
transitions  are  those  referring  to  a  Amf=  +1.  As  the  B-field  is  increased  or 
decreased  the  frequency  difference  between  the  transitional  frequencies  will 
increase  or  decrease,  respectively.  The  association  of  the  value  of  transitional 

frequencies  will  be  discussed  in  more  detail  in  the  next  section. 

An  electromagnetic  wave  applied  to  the  absorption  cell,  set  at  the 

identical  frequency  of  the  desired  transition,  will  stimulate  the  electron  to  make 
the  transition,  whereby,  the  resonance  frequency  (in  this  case  6.835  GHz)  is 
produced.  This  is  the  fundamental  principle  behind  the  operation  of  the  atomic 

clock. 
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Rubidium  Hyperfine  Energy  Levels 


The  theoretical  energy  level  equation  associated  with  the  2Si/2i  FI  and  F2 
hyperfine  splitting  can  be  calculated  from  the  following  equation  (2:29): 

E(F,  mf)=  -[(AW/4)  (21+1)]  -  [gi  (iB  B0  mf] 

+  1/2  [(AW/4)  (21+1)]  [1  +  (4mf/2I+l)K  +  K2]1/2  2-34 


where 

AW  =  4/3  [  (m  c2  a4  g;)/n3)]  [  1  +  (m/  Mp)  ]  '3  2-35 

and 

K  =  [  (gj  +  gi)  M-B  Bg]/AW(2I  +  1)  /  4]  2-36 

Most  of  the  constants  for  rubidium  have  been  discussed  in  the  previous  sections, 
but  will  be  repeated  here  for  clarity  (2:25): 

m  (electron  mass)  =  9.10910  x  10'28  g 
Mp  (proton  mass)  =  1.67252  x  10* 24  g 
g;  (Lande  g  factor  for  nuclear  spin)  =  .000995141 

gj  (Spectroscopic  splitting  factor)  =  2.002331 

mB  (Bohr  magneton)  =  9.2732  x  10'21  erg/gauss 

c  (speed  of  light)  =  2.9917925  x  108  m/sec 
a  (fine  structure  splitting  constant)  =  .00729720 

I  (nuclear  spin  magnetic  moment)  =  3/2 

B0  =  the  applied  magnetic  field  (variable) 
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With  the  exception  of  the  variable  B0,  the  only  other  factor  in  the  energy 
equation  that  is  not  a  constant  is  ntf.  As  was  seen  in  Figure  2.12,  the  F2  level  has 

five  mf  components,  whereas,  the  FI  level  has  only  three.  The  plus  or  minus 

factor  in  Equation  2.34  relates  to  the  energy  value  calculations  for  the  F2  and  F] 
hyperfine  lines,  respectively.  A  plot  of  the  theoretical  hyperfine  energy  states 
for  87Rb  is  represented  in  Figure  2.13  for  magnetic  field  values  between  0  and  1 
Gauss.  These  F2  and  FI  hyperfine  energy  values  have  also  been  provided  in 
Table  2.1  and  Table  2.2,  respectively. 

Note  that  on  Figure  2.13  an  energy  bias  was  removed  from  both  FI  and  F2 
for  plotting  purposes,  since  the  actual  energy  differentials  are  in  the 
neighborhood  of  10'21  ergs.  For  FI,  subtract  2.83  x  10' 17  ergs  from  plotted 

points  for  actual  energy  value,  and  for  F2,  add  1.7  x  10' 17  ergs. 

From  the  energy  values  listed  in  Tables  2.1  and  2.2,  the  theoretical  value 
for  the  center  frequency  (u0)  can  be  evaluated  using  the  F2;  mf=0  to  FI;  mf=0 

transitions  and  applying  Equation  2-8.  Therefore,  the  theoretical  value  for  87Rb 
(called  "Ohfs)  *s  calculated  at  6.83468261(43)  GHz  (2:24).  The  parenthesis  indicates 

the  uncertainty  digits  based  on  the  number  of  significant  digits  available  for 
the  constants. 

As  Figure  2.13  indicates,  the  energy  values  for  both  the  F2  and  FI  states  at 
the  lower  magnetic  fields  appear  to  be  linear.  Again,  the  non-linearity  occurs 
at  much  higher  B  field  values.  This  non-linearity  can  be  seen  in  Figure  2.14,  as 
the  magnetic  field  is  increased  to  4500  Gauss. 

Since  the  investigation  in  this  study  will  only  be  concerned  with  the 
linear  regions  of  the  hyperfine  splitting,  there  exists  a  linear  relationship 
between  all  the  Zeeman  transitions.  Historically,  rubidium  clock  frequency 
analysis  uses  just  one  of  the  six  possible  Zeeman  transitions  (excluding  the 
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S8o 


Table  2.1  87Rb  F2  Hyperfine  Energy  Values 
for  Coil  Magnetic  Field  Levels  From  0  to  1  Gauss 


Gauss 


B-FIELD 


mf=+2 


1 .69819159 


1 .69823794 


1 .69828429 


10' 17  Ergs 


mf=+1  mf=0 


1.69816842 


1.69819159 


1.69821477 


1.698377 


1 .69842335 


1.6984697 


1.69851605 


1.6985624 


1.69860875 


1.6986551 


1.69870145 


1.6987478 


1.69879-115 


1.6988405 


1.69888686 


1.69893321 


1.69897956 


1 .69902591 


1.69907226 | 


1.69826113 


1.69828431 


1.69830749 


1.69833067 


1.69814524 


1.69814524 


1.69814525 


1.69814525 


1.69814525 


1.69814526 


1.69814527 


1.69814527 


1.69837703 


1.69840021 


1.6984234 


1.69844658 


1.69846977 


1.69849295 


1.69851614 


1.69853933 


1.69856252 


1.6985857 


1.69860889  I 


1.69814529 


1.6981453 


1.69814531 


1.69814532 


1.69814534 


1 .6981 4535 


1.69814536 


1.69814538 


1.6981454 


1.69814541 


1.69814543 


1.69812207 


1.69809889 


1.69807572 


1.69805255 


1.69802937 


1.6980062 


1.69798303 


1.69795986 


1.69793669 


1.69791352 


1.69789036 


1.69786719 


1.69784402 


1.69782086 


1.69779769 


1.69777453 


1.69775136 


1.6977282 


1.69770504 


1.69768188 


1.69809889 


1.69805254 


1.69800619 


1.69795984 


1.69791349 


1.69786714 


1.69782079 


1.69777444 


1.69772808 


1.69768173 


1.69763538 


1.69758903 


1.69754268 


1.69749633 


1.69744998 


1.69740363 


1.69735728 


1.69731093 


1.69726458 


1.69721823 


Table  2.2  87Rb  FI  Hyperfine  Energy  Values 
for  Coil  Magnetic  Field  Levels  From  0  to  1  Gauss 


Gauss 

10-17  Ergs 

B-FIELD 

mf=-1 

!  mf=0 

mf=+1 

-2.8302188 

-2.830242071-2.83026534 

-2.83019554 

-2.83024207  -2.83028861 

-2.83017227 

-2.83024207  -2.83031188 

-2.830149 

-2.8302420 8  -2.83033515 

-2.83012574 

-2.83024208  -2.83035842 

-2.83010248 

-2.83024209  -2.83038169 

-2.83007921 

-2.83024209  -2.83040496 

-2.83005595 

-2.8302421  -2.83042823 

-2.83003269 

-2.83024211  -2.83045151 

-2.83000943 

-2.83024212  -2.83047478 

-2.82998617 

-2.83024213  -2.83049806 

-2.82996291 

-2.83024214  -2.83052133 

-2.82993965 

•2.83024215  -2.83054461 

-2.82991639 

-2.83024216  -2.83056789 

-2.82989313 

-2.83024218  -2.830591 17 

-2.82986988 

-2.83024219  -2.83061444 

-2.82984662 

-2.83024221  -2.83063772 

-2.82982337 

-2.83024222  -2.830661 

-2.82980011 

-2.83024224  -2.83068429 

-2.82977686 

-2.830242261-2.83070757 
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Figure  2.14  Plot  of  Theoretical  87Rb  F2  and  FI  HyperFine 
Energy  Values  for  Large  Magnetic  Fields 
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center  frequency),  and  is  referred  to  as  the  primary  Zeeman.  This  primary 
Zeeman  frequency  (t)z),  is  determined  from  the  the  F2;  mf=-l  to  FI;  mf=-l 
(Amf=0)  transition.  The  actual  value  of  the  frequency  is  not  used,  but  is 
represented  as  the  difference  from  the  primary  center  resonance  frequency  1)0. 
Theoretically,  there  is  a  linear  relationship  between  these  two  frequencies  (at 
low  B  field  settings)  that  can  be  calculated.  The  only  variable  that  exists  now 
from  Equation  2-34  is  the  magnetic  field  value  (B0),  since  the  transition  level 
(mf  =  -1)  for  this  primary  Zeeman  frequency  (\)z)  has  now  been  defined. 

A  similar  theoretical  relationship  can  be  made  for  the  determination  of 
the  center  resonance  frequency  (Do)  which  also  depends  on  the  value  of  the 
magnetic  field.  As  UQ  is  determined  by  the  F2;  mf=0  to  FI;  mf=0  hyperfine 
transitions,  then  a  closer  look  at  Equation  2-34  will  allow  for  certain  terms  to  go 
to  zero,  since  mf=0.  A  remaining  term  will  be  the  K  term  established  in  Equation 
2-36.  Since  B0  is  contained  in  the  evaluation  of  the  term  K,  and  the  remainder  of 
the  terms  are  constants,  then  \)0  can  also  be  represented  in  a  more  simplified 
manner  as  can  the  primary  Zeeman  frequency.  This  leads  to  a  simplified 
theoretical  expression  for  both  the  center  and  primary  Zeeman  frequencies: 

U0  =  UHFS  +  K0  B02  2-37 


and 


VZ  =  KZB0  2-38 

where  K0  is  1401940  Hz/Gauss  and  Kz  is  575.14  Hz/Gauss. 
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Both  of  these  observations  can  be  seen  in  Figures  2.15  and  2.16,  which 
plots  \)0  and  \)z,  respectively,  relative  to  varying  the  magnetic  field  between  0 

and  1  Gauss.  The  data  for  these  plots  are  represented  in  Table  2.3,  which  applies 
the  energy  values  from  each  of  the  hyperfine  lines  indicated  in  Table  2.1  and 
Table  2.2  to  the  frequency  calculation  from  Equation  2-8.  The  frequency 
associated  with  each  allowable  energy  transition  is  supplied,  but  the  individual 
transitions  of  interest  are  those  for  the  center  frequency  D0  (F2;  mf=0  to  FI;  mf=0 

transition)  and  the  frequency  (F2=  mf=-l  to  FI;  mf=-l  transition)  which  is  used 
to  determine  the  primary  Zeeman  frequency  Uz.  The  value  of  uz  is  also  provided 

in  Table  2.3,  which  has  incorporated  the  frequency  difference  calculation  from 

D0  as  discussed  above. 

Helmholtz  Coil  and  B-Field  Theory 

The  absorption  cell  B-field  used  in  the  RFS-10  design  is  a  Helmholtz  coil. 

A  representative  view  of  the  coil  configuration  used  in  the  RFS-10  rubidium 

clock  can  be  seen  in  Figure  2.17.  This  configuration  consists  basically  of  two 
circular  electrical  coils  of  the  same  radius,  a  common  axis,  and  separated  by 
some  arbitrary  distance  (depending  on  desired  results).  The  use  of  the 
Helmholtz  has  played  an  important  role  in  scientific  design  since  it  can  produce 

a  relatively  uniform  magnetic  field  over  a  small  region  of  space.  The 
determination  of  the  B-field  can  be  accomplished  by  knowing  the  dimensions  of 
the  coil,  the  number  of  turns  (N)  in  each  of  the  circular  coils,  and  the  current 
(1)  applied  to  each  of  the  coils. 

The  magnetic  induction  at  any  point  along  the  Z  axis  (boresight)  of  the 
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Figure  2.14  87Rb  Center  Frequency  Value  as  a  Function 
of  Varying  Magnetic  Field  From  0  to  1  Gauss 
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Figure  2.15  87Rb  Primary  Zeeman  Frequency  Value  as  a  Function 
of  Varying  Magnetic  Field  From  0  to  1  Gauss 
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Table  2.3  Theoretical  Fiequency  Values  for  All  Allowable  Hyperfine 
Energy  Transitions  for  ^Rb  jn  a  Magnetic  Field  from  0  to  1  Gauss 
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coil  can  be  evaluated  from  the  following  equation  (16:154): 

Bz(z)  =  (N  H0  I  a2/2)  {  [  1/  (z2  +  a2)3/2  ]  +  [  1  /  ((b-z)2  +  a2)3/2]  }  2-39 


where 


N  =  42  turns 

Ho  =  1.2566371  x  10’6  H/m 
z  =  1.3208  cm 
a  =  1.53671  cm 
b  =  .381  cm 
I  =  variable  current 

A  theoretical  value  for  the  B  field  within  the  Helmholtz  coil  can  be  seen 
graphically  in  Figure  2.18.  The  current  value  (I)  for  this  plot  was  selected  at  4.5 
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milli-amps.  The  average  B  field  within  the  absorption  cell  was  calculated  using 
a  least  squares  fit  to  the  actual  data  (4:496). 

There  are  still  some  non-uniform  aspects  to  the  Helmholtz  coil.  A  least 
squares  fit  was  applied  to  the  curve  to  determine  the  average  value  of  B0 
throughout  the  absorption  cell.  This  average  is  the  value  of  B0  during 
applications  of  Equations  2-31  through  2-34.  The  value  of  B0  was  plotted  in 
Figure  2.17  using  only  one  current  (I)  value,  but  the  shape  of  this  curve  is 
consistent  with  all  others. 


Figure  2.17  Calculated  Magnetic  Field  (B0)  for  the  Rubidium  Model  RFS-10 
Internal  Helmholtz  Coil  at  a  Current  of  4.5  Milliamps 
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III.  Experimental  Procedure 


Measurement  System 

The  C-field  experiment  was  performed  in  a  laboratory  on  the 
commercial  double-cell  Rb  frequency  standard,  EG&G  Model  RFS-10.  This 
standard  was  modified  to  allow  access  to  the  C-field  coil  wires  and  the  microwave 
power  source.  Figure  3.1  is  a  block  diagram  of  the  measurement  system.  Both  of 
the  parameters  that  are  varied,  namely  the  C-field  current  and  the  microwave 
power,  are  computer  controlled;  the  current  is  set  by  a  precision  constant- 
current  generator,  while  the  microwave  power  is  changed  by  using  an 
electronic  switch  to  change  the  resistance  of  a  bias  resistor  on  the  step- 
recovery  diode.  The  bias  resistors  were  chosen  to  change  the  microwave  power 

by  +1.3,  -1.1,  -2.2,  and  -4.0  dB  with  reference  to  the  nominal,  factoiy-set  power 
level  P0.  The  6.835-GHz  microwave  power  to  the  cell  was  sampled  by  coupling 

with  a  two-turn  coil  on  the  6.835-GHz  coaxial  signal  line.  The  bias  resistors  were 
chosen  by  using  a  Hewlett-Packard  (H-P)  model  8566  spectrum  analyzer  to 
observe  the  change  in  microwave  power  on  this  coupled  signal. 

The  entire  measurement  system  is  controlled  by  an  H-P  series  300 
computer,  which  also  acquires  and  processes  the  data.  Figure  3.2  is  a  block 
diagram  of  the  frequency  measurement  system.  The  frequency  reference  for 
both  the  Fluke  synthesizer  and  the  H-P  counter  is  an  H-P  model  5061A-004 
Cesium  (Cs)  frequency  standard. 

This  procedure  was  accomplished  following  calibration  and 
specification  steps  identified  by  the  manufacturer.  Steps  were  carefully 
taken  to  ensure  no  damage  would  occur  to  the  test  article  and  the  test  set-up 
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equipment.  After  confident  in  the  proper  operation  of  all  the  equipment, 
the  following  steps  were  taken  in  sequence  to  conduct  the  experiment. 

1.  Set  the  C- field  current  at  some  low  value  (typically  2  to  6  mA) 
and  the  microwave  power  at  some  value  (e.g.  at  the  nominal  value  PG). 

2.  Measure  the  beat  frequency  over  some  long  averaging  time  T 
(typically  500-1000  sec). 

3.  Change  the  microwave  power  level  [e.g.  to  (P0  -  1.1  dB)). 

4.  Measure  the  beat  frequency  over  T  again. 

5.  Increase  the  C-ficld  current  by  some  predetermined  amount 
(typically  0.5  mA). 

6.  Measure  the  beat  frequency  over  T  again. 

7.  Change  the  microwave  power  back  to  the  initial  value. 

8.  Repeat  steps  2  through  7  until  the  final  C-field  current 
(typically  14  to  20  mA)  is  reached. 

To  determine  the  functional  relationship  between  the  C-field 
current  and  the  Zeeman  frequency,  the  microwave  frequency  is  swept  over 
approximately  1  MHz,  centered  about  the  main  Rb  resonance  line.  The  output 
of  the  standard’s  dc-coupled  current-to-voltage  converter  is  then  plotted  as  a 
function  of  frequency.  Figure  3.3,  for  a  C-field  current  of  4.5  milli-amps,  is  a 
typical  plot.  This  plot  shows  the  main  Rb  transition  state,  as  well  as  the  four 
O  transitions  and  the  two  K  transitions  (the  71  transitions  are  used  to  define 
the  Zeeman  frequency  \)z).  It  should  be  noted  that  the  0-to-39.53-mV 
ordinate  in  Figure  3.3  rides  on  top  of  a  5  Volt  bias,  which  is  filtered  out  with  a 
precision  low-noise  floating  voltage  source.  The  Zeeman  frequency  is  read 
from  the  plot  in  Figure  3.3.  The  measurement  is  then  repeated  for  different 
B-field  currents. 
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Figure  3.3  Zeeman  Frequency  Plot  for  RFS-10 
Rubidium  Clock  at  Nominal  Factory  Setting 
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Measurement  Results 


Figure  3.4  shows  the  summary  of  results  of  measurements  made 
on  the  commercial  Rb  standard  for  changes  in  the  microwave  power  level  of 
+  1.3,  -1.1,  -2.2,  and  -4.0  dB.  Each  data  point,  which  represents  the  difference 
between  two  1000-sec  averages,  is  calculated  as  the  difference  in  output 
frequency  between  the  frequency  at  the  higher  power  and  that  at  the  lower 
power,  both  powers  being  normalized  to  the  nominal  output.  In  other  words, 

y  =  (uh  -  Di)  /  10  Mhz  3-1 

where  \)h  is  the  average  output  frequency  for  the  higher  microwave  power 
and  Uj  is  the  average  output  frequency  for  the  lower  microwave  power. 

From  Figure  3.4  we  see  that  the  maximum  frequency  change  for  the  +1.3  and 
-1.1  dB  data  is  about  3  x  lO'1*,  that  for  the  -2.2  dB  data  is  about  6  x  10"^,  and 
that  for  the  -4.0  dB  data  is  about  1.4  x  10'^.  As  a  function  of  the  change  in 
microwave  power,  the  maximum  frequency  change  is  about  2.6  x  10'^/dB. 

The  three  curves  have  a  similar  shape,  namely  one  that  is  fairly 
flat  for  Zeeman  frequencies  between  100  and  300  kHz,  then  decreasing 
monotonically  for  Zeeman  frequencies  between  300  and  850  kHz.  In  separate 
tests,  data  was  obtained  for  the  region  of  low  Zeeman  frequency.  Figure  3.5 
shows  the  results  of  those  tests  for  the  same  microwave  power  changes  as  in 
Figure  3.4.  From  the  curves  in  Figure  3.5,  it  is  seen  that  there  are  zero 
crossings  (indicating  no  frequency  shift)  at  about  13  and  37  kHz.  The 
operation  of  the  clock  at  one  of  these  points  would  result  in  zero  sensitivity 
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Figure  3.5  Fractional  Frequency  Changes  as  a  Function 
of  Small  B-Field  for  Microwave  Changes  of  +1.3,  -1.1,  -2.2,  and  -4.0  dB 
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of  the  clock  frequency  to  microwave  power  changes  and  may  result  in 
improved  long-term  clock  stability,  as  had  been  observed  by  De  Marchi  in 
his  investigations  of  cesium  frequency  standards  (12:57).  The 
manufacturer's  B-field  current  setting  for  this  clock  was  4.5  mA,  which  gave 
the  primary  Zeeman  frequency  of  approximately  191  kHz.  Operating  the 
clock  at  lower  B-fields  would  also  have  the  advantage  of  decreasing  clock 
sensitivity  to  B-field  current;  i.e.,  a  change  in  B-field  current  results  in  a 
smaller  change  in  output  frequency  when  the  B-field  current  is  small. 

A  plot  of  the  rubidium  clock  center  frequency  span  measured  at 
the  nominal  factory  magnetic  field  setting  (4.5  mA)  and  nominal  microwave 
power  setting  (PQ)  is  represented  in  Figure  3.6.  The  center  frequency  value 
(U0)  in  addition  to  the  frequency  span  at  the  half  power  point  is  indicated  on 
the  plot.  The  cell  lamp  output  voltage  indicated  on  the  y  axis  of  the  plot 
refers  to  a  relative  voltage  change  from  the  nominal  operating  output 
signal.  This  data  rides  on  1.2  Volt  cell  lamp  output  voltage,  which  has  been 
subtracted  for  the  purposes  of  plotting  these  relatively  small  (milli-volt) 
signal  changes.  The  same  data  was  collected  for  different  microwave  power 
points  which  include:  P0-  ldB,  P0-  3dB,  P0-  4dB,  P0+  ldB,  P0+  2dB,  and  P0+  3dB. 

The  results  of  these  measurements  are  indicated  in  Table  3.1.  Using  the  data 
from  this  table,  plot  of  the  half-power  linewidth  frequency  as  a  function  of 
microwave  cavity  power  measured  from  the  nominal  P0  is  represented  in 

Figure  3.7.  Also,  Figure  3.8  plots  the  relative  difference  in  center  frequency 
(U0)  measured  against  the  same  values  of  microwave  power.  There  is  an 

obvious  large  shift  in  both  half-power  frequency  width  and  center 
frequency  when  the  power  is  reduced  below  -  3dB.  Between  +  3dB,  the  half- 
power  linewidth  frequency  remains  fairly  constant  at  about  2  kHz,  and  the 
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Figure  3.6  Center  Frequency  (V0)  Line  Width  Plot  for  RFS-10  Rubidium 
Clock  Measured  at  Nominal  Factory  Power  Setting  (P0) 
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Table  3.1  Center  Frequency  (U0)  and  Half-Power  Linewidth  Frequency 
Measurements  for  Different  Microwave  Power  Values 


CAVITY  HALF  POWER  HALF-POWER  CENTER 

POWER _ SIGNAL  STRENGTH  FREQUENCY  WIDTH  _  FREQUENCY 


+  3  dB 

2.385  mV 

2112.(5)  Hz 

6834686187.(9)  Hz 

+  2dB 

2.575  mV 

2191.(1)  Hz 

6834686141.(3)  Hz 

+  1  dB 

2.445  mV 

2158.(6)  Hz 

6834686171.(6)  Hz 

Po 

2.350  mV 

2114.(2)  Hz 

6834686094.(8)  Hz 

-  1  dB 

2.185  mV 

1915.(5)  Hz 

6834686187.(1)  Hz 

-  3  dB 

1.640  mV 

1580.(4)  Hz 

6834686186.(2)  Hz 

-  4  dB 

0.215  mV 

6261. (3)Hz 

6834686020.(5)  Hz 

Figure  3.7  Measured  Rubidium  RFS-10  Half-Power  Linewidth  Frequency 
as  a  Function  of  Changing  Microwave  Power 
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Figure  3.8  Measured  Rubidium  RFS-10  Center  Frequency  Differences 
as  a  Function  of  Changing  Microwave  Power 


center  frequency  shift  is  always  in  the  positive  direction  at  about  90  Hz. 

Sweeping  the  microwave  power  from  lowest  to  highest  possible 
values  allowed  on  the  clock  adjustment,  measurements  of  the  lamp  output 
signal  voltages  were  taken  and  plotted  in  Figure  3.9.  The  selected  power 
points  discussed  above  are  annotated  on  the  plot  to  provide  for  a  point  of 
reference.  The  y-axis  of  the  pic'  provides  for  a  relative  vice  absolute  change 
in  lamp  signal  voltage  as  compared  to  the  value  observed  from  the  nominal 
factory  setting. 
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IV.  Analysis  of  Results 


Hyperfine  Theory  Versus  Experimental  Results 

Evaluation  of  the  theoretical  model  of  the  87Rb  hyperfine  transition 
frequencies  compared  to  the  measured  values  obtained  during  this  study 
requires  an  accurate  measurement  or  calculation  of  the  clock’s  coil  magnetic 
field.  As  discussed  in  Chapter  2  (pages  2-27  through  2-33),  the  theoretical  model 

for  determining  both  the  rubidium  center  frequency  (\)0)  and  primary  Zeeman 
frequency  (Uz)  rely  on  the  knowledge  of  the  magnetic  field  (B0).  Since 
measurement  data  can  only  be  obtained  from  recording  the  magnetic  field 
current,  then  alternate  methods  are  required  to  determine  an  accurate  magnetic 
field  value.  Using  the  Helmholtz  coil  model  to  calculate  the  magnetic  field  would 
not  be  appropriate  since  this  model  has  not  yet  been  validated  through  a 
theoretical  comparison.  However,  determining  the  magnetic  field  through  the 
use  of  Equations  2-34  through  2-38  will  establish  a  theoretical  baseline  to 
compare  measurement  results  against. 

One  set  of  measured  values  for  magnetic  field  current,  center  frequency 
and  primary  Zeeman  frequency  are  provided  below  which  were  used  during 
this  analysis,  however,  other  measured  values  could  have  also  been  used  to 
obtain  the  same  results. 

Coil  current,  1  =  4.50(0)  mA  (+  0.005) 

Zeeman  Frequency,  vz=  191672.(0)  Hz  (+0.1  Hz) 

Center  Frequency,  U0  =  6834686094.(8)  Hz  (+0.1  Hz) 
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The  above  Zeeman  frequency  is  applied  to  Equation  2-38  resulting  in: 

B0(t)  =  0.136719(1)  Gauss 

(where  (t)  indicates  theoretical).  The  last  digit  in  parentheses  indicates  the 

precision  uncertainty  term  based  on  the  accuracy  of  the  Zeeman  frequency 
used  (+  0.1  Hz). 

Applying  the  center  frequency  measurement  identified  above  to  Equation 

2-37  provides  the  following  magnetic  field  calculation  (DhfS  =  6.83468261(43) 

GHz  -  reference  page  2.28): 

B0(t)  =  0.049836(7)  Gauss 

There  is  approximately  a  factor  of  3  difference  between  these  two 
calculations.  Since  the  first  calculation  of  the  magnetic  field  above  is 
independent  of  the  center  frequency  (from  Equation  2-38),  then  using  the 
theoretical  energy  equation  (2-34)  should  help  in  determining  the  most 
accurate  value  for  the  magnetic  field.  Referencing  the  data  in  Tables  2-1 
through  2.3  and  using  the  value  of  B0  =  0.136719(1)  Gauss  from  above,  a 

theoretical  value  for  the  primary  Zeeman  frequency  can  be  obtained.  This 

Zeeman  is  determined  through  the  F2;  mf=-l  to  FI;  mf=-l  energy  transition.  The 
results  are  as  follows: 

\)z(t)  =  191673.(16)  Hz 
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There  is  approximately  a  1  Hz  difference  between  this  calculated  value  and  that 
which  was  measured.  This  good  comparison  indicates  that  this  value  of  the 
magnetic  field  is  accurate  to  determine  the  validation  of  the  Helmholtz  coil 
model. 

Also  using  the  theoretical  energy  equation  (2-34)  and  Tables  2-1  through 
2-3,  a  value  for  the  center  frequency  can  be  calculated  using  the  F2;  mf  =0  to  FI; 
mf=0  transition.  Again  the  same  value  (0.139719(1)  Gauss)  will  be  used.  The 
resultant  center  frequency  is 

U0(t)  =  6834682624.(0)  Hz 

The  difference  between  this  value  and  the  frequency  measured  is  3470.8  Hz.  The 
primary  factors  related  to  this  center  frequency  shift  are  related  to  both 

rubidium  clock  lamp  and  the  absorption  cell  temperatures  (6:507).  Since  the 

measurement  system  used  did  not  monitor  these  data  points  through  the 
experiment,  a  comparison  of  these  results  to  theoretical  models  could  not  be 
made.  Applying  the  known  magnetic  field  to  Equation  2-34  does  not  provide  any 
useful  information.  Center  frequency  shifts  associated  with  low  magnetic  fields 
(100  to  200  mill-Gauss)  are  on  the  order  of  200  to  300  Hz.  However,  even  though 
the  center  frequency  may  be  shifted,  the  primary  Zeeman  frequency 
relationship  is  very  stable  and  is  not  affected  by  the  same  environmental 
changes.  The  magnetic  field  value  determined  from  the  Zeeman  frequencies  are 

commonly  used  today  in  highly  accurate  space  magnetic  field  experiments. 
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Magnetic  Field  Model 


The  approach  in  evaluating  the  Helmholtz  model  also  uses  Equation  2-38. 
As  discussed  in  the  last  section,  accurate  absorption  cell  current  and  Zeeman 
frequency  measurements  were  obtained  through  the  experiment.  Applying  the 
current  values  (I)  measured  to  Equation  2-39  would  result  in  a  value  for  the 
magnetic  field  (B0(m».  Substituting  this  value  and  the  corresponding  measured 

Zeeman  frequency  into  Equation  2-38,  would  result  in  a  determination  of  the 
constant  term  Kz(m).  This  value  can  then  be  compared  to  the  theoretical 
Zeeman  frequency  constant,  Kz(a).  However,  as  seen  in  Figure  2.17  (page  2.37), 
there  is  no  constant  value  obtained  in  using  the  Helmholtz  model.  As  would  be 
expected,  the  magnetic  field  is  weaker  near  the  ends  of  the  absorption  cell  and 
strongest  in  the  center.  A  least  squares  fit  was  applied  to  each  curve  associated 
with  Figure  2.17,  so  that  a  constant  value  of  B0  could  be  used  for  each  Zeeman 

frequency  measured.  Dividing  the  measured  Zeeman  frequency  by  this 

constant  (B0)  value  will  determine  Kz(m).  Zeeman  frequency  measurements 

were  obtained  from  current  values  varying  from  3  to  20  milli-amps.  The  results 
of  this  comparison  can  be  seen  in  Figure  4.1. 

The  results  seem  to  indicate  a  fairly  large  difference  (-16274  Hz/Gauss) 
between  the  two  Kz  terms.  However,  the  error  associated  with  the  current 
measurement  will  propagate  through  the  equations  into  an  error  associated 
with  the  value  of  Kz.  The  accuracy  of  the  current  measurement  system  used  is  +. 
0.005  milli-amps.  Translating  this  accuracy  into  the  calculations  which  made  up 
Kz(m)  provides  an  error  bound  of  +16967  Hz/Gauss.  Therefore,  the  measured 
results,  including  the  error,  provide  confidence  that  Kz(m)  is  within  acceptable 
limits  of  the  theoretical  value.  This  error  bound  was  also  determined  for  the 
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Figure  4.1  Measured  Versus  Theoretical  Comparison 
of  Zeeman  Frequency  Constant  (Kz) 


accuracy  of  the  measured  BQ  term  and  calculated  as  +  0.15  milli-Gaus.  This 
accuracy  is  sufficient  when  dealing  with  magnetic  field  values  in  the  hundreds 
of  milli-gauss,  but  not  sufficient  to  evaluate  the  subtle  effects  of  the  measured 
hyperfine  states  at  very  low  (<10  milli-gauss)  ranges. 

As  shown  in  Figure  3.5,  there  exists  two  magnetic  field  values  (9  and  27 
milli-Gauss  (+  0.15),  which  resulted  in  minimizing  output  frequency  shifts.  This 
led  to  investigating  the  effects  on  the  Zeeman  frequencies  at  these  low  magnetic 
field  ranges  to  identify  potential  reasons  for  this  to  occur.  The  magnetic  field 
values  where  the  zero  crossings  occur  relate  to  current  settings  in  the  range  of 
0.2  to  0.9  milli-amps. 

To  investigate  this  area,  the  same  approach  as  before  was  taken  to 
determine  the  Zeeman  constant  Kz(m),  except  at  lower  coil  current  values 
(0  to  3  milli-amps).  Figure  4.2  represents  the  results  of  calculating  Kz(m)  at 
these  lower  current  values.  A  .01  milli-amp  increment  was  used  to  plot  the  data. 
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The  Kz(a)  value  is  not  plotted  on  this  chart  since  the  scale  would  distort  the  data 
by  compressing  the  y  axis  scale,  and  masking  some  of  the  non-linear  data.. 

These  results  seem  to  indicate  indicate  that  there  is  a  non-linear 
relationship  between  the  Zeeman  frequency  and  the  resonance  center 
frequency  at  both  large  magnetic  fields  (>1000  Gauss),  and  very  small  magnetic 
fields  (<100  milli-Gauss).  Theoretically,  the  Zeeman  frequency  determination  is 
expected  to  be  linear,  but  by  some  effect  of  the  cell  temperature  or  light 
temperature,  a  new  component  is  factored  into  the  magnetic  field  relationship. 

A  mathematical  fit  was  applied  to  this  non-linear  region  which  provided 
fairly  accurate  results  as  compared  to  the  measured  data.  The  error  per  data 
point  as  well  as  average  error  for  all  the  points  are  provided  in  Table  4.1.  The 
largest  error  to  the  fit  existed  at  the  current  value  of  0.1  milli-amps.  In  addition, 
since  there  is  a  +  0.005  milli-amp  uncertainty  in  the  measured  value,  there  will 
be  a  small  additional  error  in  the  difference  calculation.  As  can  be  seen  in  Table 
4.1,  the  measured  data  matches  the  mathematical  curve  well  within  1%.  If  the 
same  results  could  be  obtained  from  other  frequency  sources  at  these  low 
magnetic  field  settings,  it  could  possibly  validate  this  non-linear  fit.  This  fit  is 
represented  by  the  following  equation: 

Kz  =  Ao  e  -  (c/c+m)  4.1 


where 

Ao  =  1385000 
c=.  0037 

m=  current  (milli-amps) 
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Figure  4.2  Measured  Zeeman  Frequency  Constant  Kz 
at  Low  Current  Values 
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This  would  indicate  that  for  practical  clock  operations  the  calculation  for 
the  Zeeman  frequency  would  be  represented  by: 

Dz  =  Ao  e  _[c/(c+m)]  Kz  bo  4.2 

Since  this  is  based  on  only  one  clock  sample,  it  would  be  premature  to 
incorporate  this  finding  for  all  rubidium  clocks.  Further  experiments  on  a 
greater  sample  of  clocks  may  validate  these  results. 


Table  4.1  Comparison  of  Zeeman  Frequency  Model 
to  Measured  Values 


Current(ma) 

e(-c/(c+m)) 

Kz 

\WKmmm 

(Kz(m)-Kz)/Kz(m) 

mmmm 

0.1 

0.9654 

1,337,076 

1,327,568 

-0.72% 

0.2 

0.9822 

1,360,398 

1,359,550 

-0.06% 

0.3 

0.9881 

1,368,451 

1,368,244 

-0.02% 

0.4 

0.9910 

1,372,533 

1,372,592 

0.00% 

0.5 

0.9928 

1,374,999 

1,375,853 

0.06% 

0.6 

0.9940 

1,376,651 

1,377,843 

0.09% 

0.7 

0.9948 

1,377,834 

1,378,647 

0.06% 

0.8 

0.9955 

1,378,724 

1,379,520 

0.06% 

0.9 

0.9960 

1,379,417 

1,380,199 

0.06% 

1 

0.9964 

1,379,972 

1,380,734 

0.06% 

1.1 

0.9967 

1,380,427 

0.08% 

1.2 

0.9970 

1,380,806 

0.07% 

1.3 

0.9972 

1,381,128 

1,382,122 

0.07% 

1.4 

0.9974 

1,381,403 

1,382,373 

0.07% 

1.5 

0.9976 

1,381,642 

1,382,590 

0.07% 

1.6 

0.9977 

1,381,851 

1,382,781 

0.07% 

AVERAGE  of 

MAGNITUDE  of 

D, -TERENCES  = 

0.00% 
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Resonance  and  Zeeman  Frequency 


My  original  hypothesis  at  the  beginning  of  this  study  was  that  there 

should  be  a  correlation  of  Zeeman  frequency  changes  with  fractional 

frequency  differences.  However,  this  was  not  the  case.  The  Zeeman 

frequencies  were  exactly  the  same  for  all  the  power  variations,  at  the  same 
magnetic  field  settings.  The  exception  to  this  rule  was  discussed  in  the  previous 
section.  Analyzing  the  Zeeman  frequency  data  at  coil  currents  from  1  to  20 
milli-amps  reflected  changes  in  the  range  of  0.04  Hz.  Since  this  is  an  area  of 
uncertainty  in  frequency  measurement,  these  results  were  not  considered  valid. 

The  other  area  to  investigate  was  any  apparent  relationship  between  the 

center  frequency  and  fractional  frequency  as  a  function  of  the  coil  magnetic 
field.  Equation  2.34  indicates  that  there  is  a  proportional  difference  in  center 
frequency  to  the  square  of  B0.  Figure  4.3  represents  a  plot  of  the  center 

frequency  as  a  function  of  magnetic  field  from  0  to  0.6  Gauss.  Figure  4.4 
represents  a  plot  of  the  measured  fractional  frequency  data  is  represented  in 
terms  of  magnetic  field,  at  a  constant  microwave  power  (P0).  The  measured 

current  values  were  converted  to  magnetic  field  values  for  purposes  of 
comparing  the  two  figures.  The  error  associated  with  the  magnetic  field,  as 
determined  earlier,  is  negligible  on  the  scale  presented  in  Figure  4.4.  A  curve 
fit  was  applied  to  Figure  4.4  using  current  (I)  as  the  variable: 

[Af/fJ  (10-9)  =  -1.6169  +  .013027  (!)  +  .075969  (I)2  4.2 


where 
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Figure  4.3  Center  Frequency  (Theoretical)  Change  as  a  Function  of  B0 


Figure  4.4  Fractional  Frequency  as  a  Function  of  B0 
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f0  =  10  Mhz  reference  output  frequency 
I  =  current  in  milli-amps 

No  correlation  was  established  between  the  expected  center  frequency 
shift  and  the  error  observed  in  the  measurements. 

Resonance  Half  Power  Linewidth 

The  differences  in  frequency  of  the  half  power  linewidth  and  center 
frequencies  were  provided  in  Table  3.1.  The  measurements  were  taken  relative 
to  the  nominal  factory  microwave  power  setting,  PQ.  The  half-power 

measurement  is  1/2  the  change  in  resonance  output  signal  voltage  measured  in 
milli-volts  (mV). 

Figure  3.9  indicates  the  obvious  decay  in  signal  output  voltage  as  a 
function  of  decreasing  cavity  microwave  power.  As  the  power  was  decreased,  it 
was  evident  through  observing  the  plots  that  the  power  broadening  effect  made 
finding  the  center  frequency  value  very  difficult.  This  can  be  seen  in  Figure 
4.5,  which  represents  the  center  frequency  plot  at  a  point  4dB  below  the 
nominal  setting.  Since  internal  frequency  adjustments  are  made  through  lamp 
intensity  changes  to  keep  the  quartz  crystal  on  the  rubidium  center  frequency, 
it  is  quite  evident  through  the  plot  that  the  light  sensor  would  need  to  be  quite 
sensitive.  The  frequency  shift  at  these  power  settings  were  the  largest,  which  is 
what  one  would  expect. 

As  the  power  is  increased  toward  P0  the  resonance  lines  become  sharper, 
and  easier  to  determine  the  center  frequency.  This  data  substantiates  the 
fractional  frequency  results  presented  in  Figure  3.4. 
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Figure  4/  Center  Frequency  (U0)  Line  Width  Plot  for  RFS-10  Rubidium 
Clock  at  4dB  Below  Nominal  Factory  Power  Setting 
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FREQUENCY  SPAN  HZ 


In  summary,  the  maximum  frequency  change  observed  for  the 
+  1.3  and  -1.1  dB  data  is  about  3  x  10"  ^  ^ ,  for  the  -2.2  dB  data  is  about  6  x  10" 1 1, 
and  that  for  the  -4.0  dB  data  is  about  1.4  x  10' m.  As  a  function  of  the 
change  in  microwave  power,  the  maximum  frequency  change  is  about  2.6  x 
10'11/dB.  The  curves  represented  in  Figure  3.4  have  a  similar  shape, 
namely  one  that  is  fairly  flat  for  Zeeman  frequencies  between  100  and  300 
kHz,  then  decreasing  for  Zeeman  frequencies  between  300  and  850  kHz. 

There  was  also  evidence  of  zero  crossings  at  Zeeman  frequencies 
of  about  13  and  37  kHz.  The  operation  of  this  clock  at  one  of  these  points 
resulted  in  zero  sensitivity  of  the  clock  frequency  to  microwave  power 
changes. 

The  manufacturer's  B-field  current  setting  for  this  clock  was 
4.5  milli-amps,  which  gave  a  Zeeman  frequency  of  about  191  kHz. 

Operating  the  clock  at  lower  B-fields  would  also  have  the  advantage  of 
decreasing  clock  sensitivity  to  B-field  current;  i.e.,  a  change  in  B-field 
current  results  in  a  smaller  change  in  output  frequency  when  the  B-field 
current  is  small. 

The  maximum  frequency  change  due  to  microwave  power 
variations  was  about  2.6  x  10' 11  per  dB  of  power  change.  At  this  B-field 
setting  an  Allan  standard  deviation  of  0.01  dB  in  the  microwave  power 
would  result  in  an  Allan  standard  deviation  in  the  frequency  of  2.6  x  10*13 
from  the  power  changes  alone. 
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V  Conclusions  and  Recommendations 


The  experimental  results  I  have  presented  here  regarding  the 
interaction  of  the  B-field  and  the  microwave  power  on  the  output 
frequency  were  conducted  on  only  one  commercial  rubidium  clock.  The 
results  have  shown  that,  in  this  particular  case,  this  clock  is  very  sensitive 
to  microwave  power  variations  or  fluctuations. 

Adjustment  of  the  coil  magnetic  field  to  the  values  observed  in 
this  experiment  to  minimize  output  frequency  changes  is  not  recommended 
at  this  time.  Further  tests  need  to  be  conducted  on  a  variety  of  clocks  to 
validate  this  data. 

The  zero  crossings  that  were  observed,  however,  may  be 
interpreted  in  terms  of  (1)  the  opposing  effect  of  the  spatially 
inhomogeneous  light  shift  and  (2)  the  C-field  gradient;  the  light  shift 
produces  a  positive  frequency  shift,  while  the  C-field  gradient  produces  a 
negative  frequency  shift  (6:506-513).  A  different  test  configuration  would 
be  required  to  validate  this  interpretation.  I  recommend  that  these  light 
shift  effects  be  investigated  experimentally  using  the  same  frequency 
clock.  Correlation  of  data  provided  here  with  those  results  could  provide 
more  insight  into  rubidium  accuracy  and  stability  for  space  and 
commercial  applications. 
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APPENDIX  A 


CALCULATION  OF  FREQUENCY  STABILITY 

Frequency  stability  in  the  time  domain  is  expressed  as  the  two-sample 
(Allan)  variance,  a y2(t). 


CTy2(t)=  l/2[yk  +  l  -yk)2 

where  the  brackets  denotes  an  average  over  a  number  of  data  points  and  yk  is 
the  average  fractional  frequency  over  time  of  the  data  point  k. 

Stability  of  a  frequency  standard  typically  exhibits  the  characteristics  of 
the  figure  below. 


Typical  Stability  Characteristics 

The  first  pan,  "I",  reflects  the  fundamental  noise  properties  of  the  standard. 
This  decreasing  behavior  continues  with  increasing  averaging  time  until  the 
"flicker  floor"  is  reached,  II,  where  the  frequency  stability  ay  is  independent  of 
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the  averaging  time.  The  third  part  of  the  curve,  III,  represents  the  stability 
deterioration  with  increasing  averaging  time.  This  typically  occurs  at  times 
ranging  from  minutes  to  days,  depending  on  the  standard. 
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